An  Investigation  of  Cavity  Resonance  and  its 
Relationship  to  Store  Force  and  Moment  Loading 


THESIS 

Christopher  J.  Coley,  Captain,  USAF 
AFIT/GAE/ENY/11-M05 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 

Wright-Patterson  Air  Force  Base,  Ohio 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  offi¬ 
cial  policy  or  position  of  the  United  States  Air  Force,  Department  of  Defense,  or  the 
United  States  Government. 


This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to  copy¬ 
right  protection  in  the  United  States. 


AFIT/GAE /ENY/1 1-M05 


An  Investigation  of  Cavity  Resonance  and  its 
Relationship  to  Store  Force  and  Moment  Loading 


THESIS 


Presented  to  the  Faculty 
Department  of  Aeronautics  and  Astronautics 
Graduate  School  of  Engineering  and  Management 
Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Gommand 
In  Partial  Fulhllment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Aeronautical  Engineering 


Ghristopher  J.  Goley,  B.S. 
Gaptain,  USAF 


March  2011 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


AFIT/GAE/ENY/11-M05 


An  Investigation  of  Cavity  Resonance  and  its 
Relationship  to  Store  Force  and  Moment  Loading 


Christopher  J.  Coley,  B.S. 
Captain,  USAF 


Approved: 


m  Col  Ricfrara  Kauffman,  PhD  (Member)  Date 


AFIT/GAE /ENY/1 1-M05 


Abstract 

A  store  which  is  released  from  an  internal  bay  is  subjected  to  a  highly  unsteady 
flowheld  which  influences  the  release  characteristics  for  the  store.  Pressure  transducer 
information  is  often  used  to  estimate  store  loading,  although  no  direct  correlation 
between  the  observed  frequencies  from  the  pressure  transducers  and  the  store  loading 
has  been  developed.  The  relationship  between  the  acoustic  modes  present  in  the  cavity 
and  the  force  and  moment  loading  on  a  store  released  from  a  bay  are  investigated 
through  a  CFD  study  using  the  OVERFLOW  2.1  solver.  The  acoustic  modes  of  the 
cavity  are  calculated  from  the  pressure  fluctuation  histories  along  the  cavity  ceiling 
and  walls.  Empty  cavity  solutions  are  compared  to  the  frequencies  predicted  by  the 
Rossiter  equation  and  to  experimental  data  from  the  WIGS  database.  Additionally, 
the  cavity  pressure  fluctuation  spectra  are  compared  to  the  force  and  moment  loading 
spectra  for  a  store  located  in  the  cavity  in  both  the  carriage  position  and  the  shear 
layer.  It  is  determined  that  the  presence  of  a  store  in  the  cavity  does  not  alter  the 
fundamental  acoustic  modes  in  the  cavity.  Also,  the  store  force  and  moment  loadings 
are  observed  to  correlate  with  the  cavity  pressure  fluctuations. 
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An  Investigation  of  Cavity  Resonance  and  its 
Relationship  to  Store  Force  and  Moment  Loading 

I.  Introduction 

Cavity  flows  have  been  the  subject  of  many  studies  since  the  1950s  [1].  The 
grazing  flow  over  cavities  is  relevant  to  a  wide  range  of  real-world  applications,  from 
gas  transport  systems  to  the  the  sunroofs  of  automobiles.  The  characteristics  of 
cavity  flow  is  an  area  of  interest  to  the  Air  Force  because  it  also  describes  the  flow 
phenomenon  experienced  in  the  region  of  an  aircraft  weapons  bay. 

Although  the  airflow  around  an  aircraft  is  unsteady,  historic  practice  has  been 
to  use  quasi-steady  techniques  wherein  the  aerodynamic  loading  on  the  store  are 
time-averaged.  These  techniques  capture  the  mean  properties  of  the  airflow  but  ne¬ 
glect  the  time-varying  fluctuations  in  the  flow.  Due  to  the  demonstrated  success  and 
repeatability  of  this  methodology  for  externally-carried  stores,  it  has  long  been  ac¬ 
cepted  that  the  release  of  an  externally-carried  store  can  be  considered  a  steady-state 
problem  [2].  This  assumption  does  not  necessarily  apply,  however,  for  the  release  of 
internally-carried  stores.  For  aircraft  with  weapon  bays,  a  shear  layer  develops  across 
the  opening  of  the  cavity.  The  interaction  of  the  shear  layer  and  its  reattachment  at 
the  aft  wall  of  the  cavity  produces  acoustic  waves  which  lead  to  highly  turbulent  flow 
in  the  cavity  and  shear  layer.  The  influence  of  the  unsteady  flow  produces  unsteady 
aerodynamic  forces  on  a  store  being  released  from  the  cavity,  possibly  affecting  the 
trajectory  of  the  store.  It  has  been  shown  that  this  unsteady  airflow  can  influence 
the  separation  of  a  store  such  that  differing  trajectories  can  result  at  different  release 
times,  even  when  all  other  flight  conditions  remain  constant  [3,4]. 

Today’s  and  tomorrow’s  military  aircraft  are  being  designed  with  an  increased 
emphasis  on  low-observability.  This  design  goal  has  led  designers  to  resort  to  internal 
carriage  of  weapons  in  order  to  reduce  radar  cross-section.  Although  internal  carriage 
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of  weapons  and  their  release  from  a  bay  is  nothing  new,  the  increasing  shift  to  smaller, 
lighter  munitions  and  ranged  weapons  has  introduced  a  new  challenge  to  the  aircraft- 
store  compatibility  engineer.  Previous  weapons  were  typically  larger  and  heavier 
and  had  high  moments  of  inertia  that  rendered  them  relatively  insensitive  to  the 
fluctuating  flow  present  in  the  region  of  a  weapon  bay.  Additionally,  it  was  common 
practice  to  use  high  ejector  forces  in  order  to  release  the  store  at  high  velocities, 
minimizing  the  store’s  time  in  the  turbulent  flow  and,  thus,  minimizing  the  impact 
on  the  trajectory.  Modern  weapons  are  more  sensitive  than  previous  weapons  to  the 
aerodynamic  loads  they  experience  as  they  are  released  from  the  aircraft  for  several 
reasons  [4]: 

•  Smaller,  lighter  stores  have  less  inertia  and  are  therefore  influenced  more  by 
fluctuating  aerodynamic  loads 

•  Modern  stores  are  designed  to  be  inherently  less  stable  in  order  to  increase 
maneuverability 

•  Stores  have  relatively  larger  control  surfaces  than  before,  making  them  more 
susceptible  to  the  unsteady  flow  environment 

•  Many  modern  stores  have  non-circular  cross-sections,  which  can  induce  rolling 
moments 

•  The  use  of  active  seeker  heads  and  guidance  systems  lead  to  attitude  constraints 
and  dynamic  limits  during  separation 

•  Sensitive  equipment  on  board  modern  weapons  oftentimes  limit  the  magnitude 
of  ejector  forces  that  can  be  used. 

Much  of  the  early  research  in  the  area  of  cavity  flows  focused  on  the  acoustic 
modes  as  it  related  to  structural  fatigue  in  the  weapon  bay.  Suppression  devices 
were  developed  as  a  result  which  served  to  quell  the  turbulent  airflow  and  hence  the 
vibrations  in  the  bay.  As  these  devices  quieted  the  airflow  in  the  bay,  they  also  served  a 
role  in  improving  store  separation  characteristics  and  predictability.  Passive  and  open- 
loop  active  suppression  devices  do  not  require  knowledge  of  the  flow  characteristics. 
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but  closed-loop  active  suppression  devices  require  the  measurement  or  estimation  of 
some  flow  quantity  [1].  As  such,  an  understanding  of  the  acoustic  response  to  cavity 
flow  is  necessary  in  order  to  design  these  devices. 

In  1970,  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  conducted  a  com¬ 
prehensive  wind  tunnel  test  program  in  order  to  characterize  the  aeroacoustic  response 
of  cavities  both  with  and  without  stores  present  [5].  From  this  study  it  was  demon¬ 
strated  that  the  presence  of  a  relatively  small  store  in  a  cavity  does  not  greatly  affect 
the  vibrational  modes  observed  in  the  cavity,  although  this  study  did  observe  that  a 
store  located  in  the  shear  layer  close  to  the  leading  edge  acted  as  a  sort  of  suppres¬ 
sion  device.  This  study  looked  only  at  the  pressure  fluctuations  in  the  cavity  and 
did  not  look  at  the  aerodynamic  forces  nor  moments  on  the  store.  From  September 
1986  through  September  1990,  Arnold  Engineering  Development  Center  (AEDC),  un¬ 
der  the  Weapons  Internal  Carriage  and  Separation  program  (WICS),  developed  an 
aerodynamic  database  for  internal  store  carriage  and  jettison  [6].  This  study  charac¬ 
terized  pressure  fluctuations  in  cavities  of  various  geometries  and  took  time-averaged 
store  forces  and  moments  for  several  stores,  but  did  not  correlate  cavity  pressure 
fluctuations  with  unsteady  loading  on  the  stores. 

Current  aircraft /store  compatibility  methodology  utilizes  time-averaged  store 
forces  and  moments  and  flowheld  data  with  a  6-degree  of  freedom  (6-DOF)  analysis 
in  order  to  simulate  expected  trajectories.  The  data  can  come  either  from  wind- 
tunnel  testing  or  from  computational  fluid  dynamics  (CFD).  The  results  from  this 
analysis  dictate  which  release  points  are  conducted  in  flight  test.  It  is  therefore  crucial 
that  the  analysis  can  accurately  predict  the  separation  characteristics  of  a  store  in 
order  to  minimize  cost  of  the  flight  test  program  and  to  ensure  aircrew  safety  during 
flight  testing.  This  requirement  drives  the  need  to  develop  a  time-accurate  predictive 
capability  for  the  analysis  of  cavity-released  stores.  In  doing  so,  it  is  important  to 
understand  the  relationship  between  the  characteristics  of  the  airflow  and  the  resultant 
trajectories.  While  it  has  been  assumed  that  cavity  pressure  fluctuations  correspond 
to  store  force  and  moment  loading  variations,  no  direct  correlation  has  been  developed. 
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If  it  can  be  determined  that  there  is  a  known  correlation  between  the  vibrational 
freqnencies  demonstrated  in  an  empty  cavity  and  the  forces  and  moments  on  a  store 
in  that  cavity,  then  a  better  predictive  capability  for  store  separation  analysis  can  be 
developed. 

The  objective  of  the  present  research  is  to  develop  a  greater  nnderstanding  of 
the  role  that  the  acoustic  modes  present  in  a  cavity  flow  have  on  the  aerodynamic 
forces  and  moments  imposed  upon  a  store  released  from  that  cavity.  A  CFD  study 
is  conducted  which  correlates  the  cavity  unsteady  pressure  data  with  store  unsteady 
force  and  moment  loading.  The  time-accurate  solutions  for  an  empty  cavity,  a  cavity 
with  a  weapon  in  the  carriage  position,  and  a  weapon  located  in  the  shear  layer  are 
investigated.  A  generic  cavity  with  a  length-to-depth  ratio  oi  L/D  =  4.5  and  a  generic 
weapon  are  utilized  in  the  study  in  order  to  allow  for  validation  with  data  from  the 
WIGS  database.  The  flight  condition  is  M  =  0.95,  which  represents  a  typical  worst- 
case  point  within  the  flight  release  envelope.  The  vibrational  modes  present  in  the 
empty  cavity  are  compared  to  those  in  the  two  cases  with  the  weapon  present  and 
with  the  aerodynamic  forces  and  moments  on  the  body  of  the  store  in  both  cases. 

The  introduction  has  laid  the  framework  of  this  research  and  the  motivation 
behind  it.  Chapter  II  provides  background  on  the  physics  of  the  cavity  flow  phe¬ 
nomenon  as  well  as  an  overview  of  previous  cavity  flow  research.  This  chapter  also 
describes  the  OVERFLOW  solver,  which  is  used  in  the  current  research.  Chapter  III 
describes  the  methodology  used  in  this  research,  detailing  the  conhgurations  which 
are  examined.  Results  are  located  in  Chapter  IV.  Lastly,  Chapter  V  is  composed  of 
the  conclusions  and  recommendations  for  future  work. 
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II.  Background 

The  present  work  seeks  to  correlate  the  pressure  measurements  found  in  cavity  flow 
to  the  unsteady  weapon  loads  imposed  upon  a  store  located  in  a  cavity.  A  general 
characterization  of  cavity  flow  is  provided  and  a  description  of  the  WIGS  wind  tunnel 
test  and  its  relevant  hndings  are  summarized.  The  OVERFLOW  2.1  solver,  which 
is  used  to  calculate  the  solutions  for  this  work,  is  described  and  its  main  features 
outlined.  Finally,  a  brief  overview  of  signal  processing,  along  with  a  few  considerations 
regarding  spectral  analysis,  are  presented. 

2.1  Characterization  of  Cavity  Flow 

The  grazing  flow  past  a  cavity  has  been  studied  with  great  interest  as  it  has  ap¬ 
plications  in  areas  including,  but  not  limited  to  automotive  components,  gas  transport 
systems,  and  aeronautical  systems.  The  underlying  issue  of  concern  is  the  occurrence 
of  high  amplitude  pressure  oscillations  as  they  lead  to  unsteady  loading  on  the  com¬ 
ponents  within  or  near  the  cavity  or  unwanted  levels  of  noise  in  the  cavity. 

A  shear  layer  is  produced  above  the  cavity  opening  when  the  approaching  bound¬ 
ary  layer  separates  at  the  leading  edge  of  the  cavity.  The  cavity  flow  is  then  char¬ 
acterized  by  the  shear  layer  behavior.  Rockwell  and  Naudascher  [7]  identihed  three 
flow  regimes  which  they  labeled  fluid-dynamic  interactions,  fluid-resonant  interac¬ 
tions,  and  fluid-elastic  interactions.  The  fluid-dynamic  regime  involves  shear-layer 
instability  amplihcation  due  to  feedback  with  interaction  of  the  shear  layer  with  the 
aft  cavity  wall.  This  phenomenon  is  depicted  in  Fig.  1.  This  interaction  is  typical  in 
low-speed  flow  past  shallow  cavities.  The  fluid-resonant  regime  couples  the  acoustic 
modes  of  the  cavity  and  the  shear  layer  over  the  cavity.  In  this  regime,  the  cavity 
oscillations  are  strongly  coupled  with  wave  effects  in  the  cavity.  This  regime  occurs 
when  the  frequencies  are  high  enough  that  the  acoustic  wavelength  is  on  the  order 
of  or  smaller  than  the  cavity  characteristic  length  [7].  This  interaction  is  typical 
in  deeper  cavities  and  cavities  located  in  high  Mach  number  flow.  The  fluid-elastic 
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regime  encompasses  flows  that  are  affected  by  the  elastic  boundaries  of  the  cavity  and 
are  typically  encountered  when  one  part  of  the  cavity  is  being  actuated. 


Boundary  Layer 


Figure  1:  Schematic  depicting  acoustic  resonance  resulting  from  grazing  cavity  flow 

Cavities  are  typically  dehned  as  shallow  or  deep  based  on  their  length-to-depth 
ratio  {L/D).  Rossiter  [8]  dehnes  shallow  cavities  as  those  having  L/D  >  4  and  deep 
cavities  as  those  having  L/D  <  4,  although  these  dehnitions  are  not  universal. 

At  low  speeds,  the  cavity  flow  can  be  described  as  being  in  the  shear-layer  mode 
or  wake  mode  [9] .  The  shear-layer  mode  occurs  when  the  shear  layer  spans  the  opening 
of  the  cavity  and  stagnates  at  the  aft  wall.  A  cavity  in  this  mode  is  also  known  as 
an  open  cavity.  In  this  classihcation,  the  flow  can  experience  both  fluid-dynamic  and 
fluid-resonant  interactions.  The  wake  mode  occurs  when  the  shear  layer  stagnates 
somewhere  prior  to  the  aft  wall  (e.g.  reattachment  with  the  cavity  base).  A  cavity 
in  this  mode  is  also  known  as  a  closed  cavity.  For  supersonic  flows  there  can  also  be 
interim  conditions  which  are  known  as  transitionally  open  and  transitionally  closed. 
The  two  different  types  of  cavity  flow  are  depicted  in  Figs.  2  and  3.  Initially  the  cut-off 
between  an  open  and  a  closed  cavity  were  thought  to  be  determined  by  the  length-to- 
depth  ratio  alone.  However,  various  sources  throughout  the  literature  cite  different 
values  of  L/D  as  the  cut-off  between  open  and  closed,  for  instance:  L/D  <  10  for 
open  and  L/D  >  13  for  closed  [10],  L/D  <  9  for  open  and  L/D  >  13  for  closed  [11], 
L/D  <  3  for  open  and  L/D  >  10  for  closed  [12].  It  is  believed  that  these  discrepancies 
are  due  to  the  freestream  Mach  number  and  boundary  layer  thickness,  in  addition  to 
L/D,  affecting  the  observed  mode  [9]. 
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Figure  2;  Depiction  of  an  open  cavity 


Figure  3:  Depiction  of  a  closed  cavity 

For  open  cavities,  the  freqnencies  of  oscillation  can  be  predicted  using  empirical 
formulae.  For  speeds  of  M  >  0.4,  the  frequencies  can  be  approximated  by  the  Rossiter 
eqnation  [8]: 


U  m  —  a 
L  4  +  M’ 

K 


(1) 


where  U  is  the  freestream  velocity,  m  is  the  mode  number,  and  Rossiter  dehned  a 
and  K  experimentally  as  0.25  and  0.57,  respectively.  Note  that  Eq.  (1)  only  considers 
longitndinal  modes  and  is  therefore  a  two-dimensional  model.  In  the  two-dimensional 
model,  it  is  assnmed  that  the  stagnating  flow  at  the  cavity  aft  wall  is  deflected  at  the 
same  frequency  as  the  vortex  shedding  at  the  leading  edge  (although  not  in  phase).  In 
the  three-dimensional  case,  however,  the  cavity  response  is  inflnenced  by  the  geometry 
of  the  cavity,  the  condition  of  the  snrfaces  of  the  cavity,  the  contents  of  the  cavity, 
and  the  characteristics  of  the  flnid  in  the  cavity.  Conpling  of  the  vortex  shedding  and 
fluid  response  is  only  approximated  by  the  two-dimensional  Rossiter  model.  While 
the  longitudinal  modes  dominate  the  flow,  it  is  plausible  that  lateral  and  vertical 
modes  also  exist  [6].  The  Rossiter  eqnation  assnmes  a  speed  of  sonnd  equal  to  the 
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freestream  speed  of  sound  and  as  a  result  produces  larger  errors  as  the  Mach  number 
increases.  Heller,  et  ah  [5]  modihed  the  Rossiter  equation  to  provide  smaller  error  as 
Mach  number  increases.  The  modihed  Rossiter  equation  is 


/ 


U  m  —  a 

1  j _ M 


(2) 


where  7  represents  the  ratio  of  specihc  heats.  Note  that  Eqs.  (1)  and  (2)  predict 
the  expected  frequencies  but  give  no  information  regarding  the  amplitude  of  these 
frequencies.  Additionally,  these  equations  apply  only  to  open  cavities,  so  it  must  be 
known  which  how  regime  is  present  for  the  cavity  in  question.  Therefore  simulation  or 
experimentation  are  required  to  completely  determine  the  characteristics  of  the  how. 


The  acoustic  frequencies  in  the  bay  can  be  related  to  the  time  history  of  the 
pressure  huctuations  P{t).  This  relationship  is  accomplished  by  computing  the  power 
spectral  density  (PSD),  which  describes  how  the  power  of  P{t)  is  distributed  with 
frequency.  The  PSD  is  computed  from  the  Fourier  transform  of  the  autocorrelation 
of  P{t).  From  the  convolution  theorem,  it  can  be  dehned  as  the  squared  modulus 
of  the  Fourier  transform  of  P{t)  [13].  It  is  common  to  report  sound  pressure  on  a 
logarithmic  scale  relative  to  the  lowest  human  audible  sound.  This  logarithmic  scale 
is  known  as  the  sound  pressure  level  (SPL).  Sound  pressure  level  is  related  to  PSD 
by  the  following  equation: 


SPL 


(3) 


where  pref  =  2  x  10  ®  Pa  is  the  value  adopted  as  the  minimum  audible  sound  pressure 
variation. 


2.2  WICS  Program 

Arnold  Engineering  Development  Center  conducted  a  wind-tunnel  test  program 
from  September  1986  through  September  1990  with  the  intention  of  developing  an 


aerodynamic  database  for  use  in  analysis  of  internal  store  carriage  and  jettison.  Dur¬ 
ing  the  course  of  this  program,  aerodynamic  forces  acting  on  a  series  of  slender  missile 
conhgurations  were  measured  with  the  missiles  mounted  in  a  spatial  grid  of  locations 
in  and  near  a  generic  cavity  in  a  flat  plate.  The  cavity  was  equipped  with  a  vari¬ 
ety  of  doors  and  acoustic  suppression  devices,  such  as  boundary-layer  spoilers  and 
aft-bulkhead  ramps.  Static  and  fluctuating  pressures  acting  on  the  plate  and  cavity 
surfaces  were  also  measured.  The  fluctuating  pressures  were  processed  into  spectra 
via  Fourier  transform  techniques,  and  the  resonant  modes  found  in  the  cavity  were 
determined.  The  experiments  were  conducted  with  a  variety  of  cavity  geometries  at 
Mach  numbers  ranging  from  0.60  to  5.04  [6].  The  data  collected  during  this  program 
has  served  as  a  source  by  which  CFD  simulations  have  been  validated.  Only  the 
experimental  setup  and  results  pertaining  to  the  collection  of  the  pressure  fluctuation 
data  for  the  empty  cavity  are  presented  here,  as  they  pertain  to  the  present  research. 

2.2.1  WICS  Experimental  Setup.  The  cavity  in  the  experiment  was  con¬ 
structed  by  building  an  opening  4  in  wide  by  18  in  long  into  a  flat  plate  16  in  wide  by 
47  in  long.  The  cavity  ceiling  could  be  installed  at  any  of  several  discrete  heights  be¬ 
tween  0  and  4  in.  Only  1.25-,  2-,  and  4-in  heights  were  used  during  the  tests,  providing 
cavities  of  length-to-depth  ratios  oi  L/D  =  14.4,  9.0,  and  4.5,  respectively.  Fluctuat¬ 
ing  pressures  were  measured  with  differential  transducers  at  up  to  45  locations:  7  on 
the  flat  plate  and  38  on  the  walls  and  ceiling  of  the  cavity.  Not  all  transducers  were 
necessarily  used  in  every  test.  Data  was  collected  at  a  rate  of  10,000  samples/sec. 
The  test  conhguration  is  illustrated  in  Fig.  4. 
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Figure  4:  Top  view  of  plate  and  cavity  setup  as  mounted  in  wind  tunnel 

Data  were  recorded  at  Mach  numbers  in  the  range  from  0.60  to  5.04.  A  nominal 
unit  Reynolds  number  of  3  x  10®//t  was  selected,  but  unit  Reynolds  number  actually 
varied  during  the  tests.  A  full  listing  of  nominal  values  of  the  flow  conditions  can  be 
found  in  the  WIGS  technical  report  [6].  The  nominal  flow  conditions  for  the  M  =  0.95 
cases  are  listed  in  Table  1.  Note  that  the  case  numbering  is  a  convention  used  here 
and  was  not  used  during  the  WIGS  program. 


Table  1:  Nominal  flow  conditions  for  WIGS  M  =  0.95  tests 


Gase 

Moo 

Pt 

psf 

Tt 

deg  R 

Roo 

ft /sec 

Qoo 

psf 

Re 

1/ft 

1 

0.95 

478 

542 

998 

169 

1.0  X  10® 

2 

0.95 

980 

545 

1000 

343 

2.0  X  10® 

3 

0.95 

1200 

550 

1008 

424 

2.5  X  10® 

4 

0.95 

1480 

551 

1008 

525 

3.0  X  10® 

2.2.2  WICS  Results.  Modal  frequency  predictions  for  the  WIGS  cavities 
were  made  using  the  modihed  Rossiter  equation  with  values  of  a  and  k  of  0.28  and 
0.56  for  the  L/D  =  4.5  and  LjD  =  9.0  cavities,  respectively.  These  values  were 
determined  during  the  WIGS  program  to  provide  superior  matching  than  the  values 
originally  posed  by  Rossiter  [6].  Gomparisons  between  the  the  predicted  and  measured 
frequencies  are  presented  in  Table  2.  The  modal  amplitudes  were  found  to  be  in  the 
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range  of  130  to  170  dB.  Note  that  the  L/D  =  14.4  cavity  is  a  closed  cavity  and 
therefore  the  modified  Rossiter  equation  does  not  apply  to  this  case.  Additionally,  it 
was  reported  that  no  modes  were  detected  during  the  experiment  for  this  cavity  [6]. 
Of  specihc  note  are  the  results  from  the  L/D  =  4.5  cavity  at  M  =  0.95,  since  this  is 
the  case  which  is  examined  in  the  present  research. 


Table  2:  Comparison  of  predicted  and  measured  modal  frequencies  from  the  WIGS 
test  program 


Mach  number 

0.60 

0.95 

1.20 

1.50 

2.00 

L/D 

Mode 

Pred. 

Meas. 

Pred. 

Meas. 

Pred. 

Meas. 

Pred. 

Meas. 

Pred. 

Meas. 

4.5 

1 

137 

137 

185 

186 

210 

225 

230 

244 

254 

303 

2 

327 

352 

441 

469 

502 

518 

549 

586 

608 

664 

3 

518 

547 

698 

732 

794 

830 

869 

938 

961 

1,064 

4 

708 

752 

954 

1,016 

1,087 

1,025 

1,188 

1,152 

1,314 

1,436 

5 

898 

967 

1,211 

1,318 

1,379 

1,338 

1,507 

1,318 

1,668 

1,797 

9.0 

1 

84 

— 

113 

— 

129 

— 

141 

— 

155 

— 

2 

274 

313 

369 

400 

421 

439 

460 

498 

509 

— 

3 

464 

— 

626 

684 

713 

762 

779 

859 

862 

625 

4 

655 

762 

883 

986 

1,005 

1,104 

1,099 

1,230 

1,216 

1,016 

5 

845 

— 

1,139 

1,270 

1,297 

1,396 

1,418 

1,592 

1,569 

1,455 

The  results  from  the  WIGS  experiment  have  been  shown  to  be  reproduceable  by 
CFD  simulation.  Nichols  [14]  investigated  the  flow  characteristics  of  the  L/D  =  4.5 
cavity  with  M  =  0.95  utilizing  Detached  Eddy  Simulation  (DES)  turbulence  mod¬ 
els.  These  solutions  were  computed  using  a  second-order  time,  third-order  upwind 
HHLEM  [15]  space  algorithm.  These  results  were  in  good  agreement  with  the  data 
for  the  time-averaged  pressure,  overall  sound  pressure  level,  and  sound  pressure  level 
spectrum.  In  general,  the  first  four  acoustic  modes  were  identifiable  in  the  predic¬ 
tions.  In  another  study,  Nichols  [16]  investigated  the  ability  of  Delayed  Detached  Eddy 
Simulation  (DDES)  turbulence  models  with  and  without  rotation  and  curvature  cor¬ 
rections  to  accurately  predict  the  flowheld  characteristics  of  both  the  L/D  =  4.5 
and  L/D  =  9.0  cavities.  In  general,  it  was  found  that  the  DDES  simulations  pro¬ 
vided  slightly  superior  results  to  those  produced  by  DES,  although  the  rotation  and 
curvature  corrections  made  no  signihcant  improvement  to  the  solutions. 
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2.3  OVERFLOW  2.1  CFD  Solver 


The  present  work  uses  the  OVERFLOW  2.1  CFD  solver.  This  code  is  a  three- 
dimensional  time-marching  implicit  Navier-Stokes  code  that  can  also  operate  in  two- 
dimensional  or  axi-symmetric  mode  which  utilizes  structured  overset  grids.  OVER¬ 
FLOW  2.1  has  several  different  inviscid  ffux  algorithms  and  implicit  solution  algo¬ 
rithms  and  includes  both  thin  layer  and  full- viscous  terms.  The  OVERFLOW  2.1 
code  includes  algebraic,  one-equation,  and  two-equation  turbulence  models,  including 
three  DES  turbulence  models  [17].  The  code  also  supports  bodies  in  relative  mo¬ 
tion  and  includes  a  6-DOF  model.  Parallelization  is  possible  through  both  MPI  and 
OpeuMP. 

2.3.1  OVERFLOW  2.1  Capabilities.  The  following  is  a  partial  listing  of  the 
capabilities  currently  included  in  OVERFLOW  2.1  [18]: 

Inviscid  ffux  algorithms: 

•  2nd/4th/6th-order  central  difference  with  smoothing 

•  Yee  Symmetric  TVD 

•  AUSM-h  upwind 

•  3rd-order  Roe  upwind 

•  3rd-order  HLLC  upwind 

•  5th-order  WENO  and  WENOM  upwind 
Implicit  solvers: 

•  ADI  Beam- Warming  block  tridiagonal  solver  with  either  central  difference  or 
upwind  Steger- Warming  ffux  Jacobians 

•  Steger-Warming  2-factor  scheme 

•  ADI  Pulliam-Chaussee  scalar  pentadiagonal  solver 

•  LU-SGS  solver 
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•  D3ADI  diagonalized  solver 

•  SSOR  solver 
Turbulence  models: 


•  Baldwin-Lomax  algebraic  model  with  wake  model 

•  Baldwin-Barth  1-equation  transport  model 

•  Spalart-Almaras  1-equation  transport  model 

•  Spalart-Almaras  DES  and  DDES  models 

•  k-ca  2-equation  transport  model 

•  SST  2-equation  transport  model 

•  SST  DES  and  SST  DDES  2-equation  transport  models 

•  SST  Multi-Scale  (MS)  2-equation  transport  model 

•  Wall  functions  are  available  for  all  of  the  transport  turbulence  models 

•  Rotation  and  curvature  corrections  are  available  for  the  1-  and  2-equation  mod¬ 
els 

•  Temperature  correction  is  available  for  the  2-equation  models. 

2.3.2  Navier- Stokes  Equation  Implicit  Solution  Procedure.  OVERFLOW 
2.1  solves  the  Navier-Stokes  in  generalized  coordinates.  The  following  treatment  is 
derived  from  the  OVERFLOW  2.1  User’s  Guide  [18].  The  Navier-Stokes  equations 
may  be  written  as 


dq  dE  OF  dG 
dt  ^  ^  dp  ^  d( 

where  q  is  the  vector  of  conserved  variables 


(4) 
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q  = 
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pv  , 

pw 
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(5) 


where  V  represents  volume.  The  linearized  Euler  implicit  form  of  Eq.  (4)  including 
sub-iterations  is  given  by 


At 


{i+e)AT  ^  i+e 


/  + 

|^(gn+l,'m 


T  T  d-qB  -\-  d(^C) 


/Aq 


f  I  j- 


g”)  -  YfjAg”  +  . 


(6) 


Here  0  =  0  for  first  order  time  differencing  and  0  =  1/2  for  second  order  time 
differencing.  An  artihcial  time  term  j  has  been  explicitly  added.  The  pseudo 

time  (Ar)  may  vary  throughout  the  flow  held  when  a  local  time  step  is  employed.  The 
artihcial  time  term  must  converge  at  each  physical  time  step  to  assure  time  accuracy. 
The  explicit  viscous  and  inviscid  huxes  are  included  in  the  term  RHS  given  by 


RHS 


dE  dF  dG 
dp  d(  ' 


(7) 


Equation  (6)  is  of  the  form  Ax  =  b,  where  the  hrst  bracketed  term  in  (6)  is  the 
left  hand  side  matrix  A  and  the  second  bracketed  term  is  the  vector  b.  This  system  of 
equations  can  be  solved  by  inversion  of  the  A  matrix.  Direct  inversion  of  the  matrix 
A  requires  considerable  computational  time  and  memory,  so  approximations  are  used 
to  expedite  the  procedure.  Equation  (6)  can  be  factored  in  space  as 


+  AAA]  [I  +  [I  +  A, = 

_  _  -|^Ag"  +  +  Error 

where  the  factorization  error  is  given  by 
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Error  =  {d^Adr^B  +  d^Ad^^C  +  dr,Bd^C) 

+  (^)^  {d^Adr,Bd^C)\  Aq^+^’^+^ 


(9) 


and  A,  B,  and  C  are  block  tridiagonal  matrices  for  structured  grids  with  central 
difference  or  first-order  spatial  upwind  implicit  flux  Jacobians.  The  A,  B,  and  C 
matrices  can  be  decomposed  into  eigenvalues  (A)  and  eigenvectors  (X)  as 

A  =  XaAaX-/ 

B  =  XbAbX^^  (10) 

C  =  XcAcXc^. 

One  method  to  solve  the  unfactored  system  of  equations  in  Eq.  (6)  is  to  use 
relaxation  procedures.  OVERFLOW  2.1  utilizes  one  such  method  with  the  Symmetric 
Successive  Over- Relaxation  scheme  (SSOR).  This  scheme  can  be  written  as 


=  (1  -  f!)  ACw  +  ^ 


-B.Aq-^El  -  BRAqftE  -  CqAqJ^E,  -  CuAq^El)  ■ 


where  the  subscripts  L  and  R  denote  the  left  and  right  blocks  of  the  tridiagonal 
matrices,  respectively.  The  overbar  indicates  a  pre-multiply  by  the  inverse  of  the 
diagonal  matrix  Ad  +  Bd  +  Cd  where  the  subscripts  D  denotes  the  diagonal  block  of 
the  tridiagonal  matrices.  The  update  level  of  Ag  is  given  by  mm.  The  scheme  uses  a 
forward  sweep  in  j  and  symmetric  sweeps  in  k  and  /.  For  a  forward  sweep  in  k  and  I, 
the  updated  levels  are  dehned  as 


mkl  =  mm  +  l,mk2  =  mm,  mil  =  mm  +  1,  m/2  =  mm.  (12) 

For  a  backward  sweep  in  k  and  /,  the  update  levels  are  dehned  as 
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mkl  =  mm,  mk2  =  mm  +  1,  mil  =  mm,  ml2  =  mm  +  1. 


(13) 


A  symmetric  sweep  consists  of  a  forward  and  a  backward  sweep.  Multiple 
symmetric  sweeps  (typically  10)  are  performed  at  each  sub-iteration.  The  relaxation 
parameter  is  typically  set  to  0.9. 

2.3.3  Overset  Grids.  OVERFLOW  2.1  only  supports  structured  grids.  One 
difficulty  that  arises  from  the  use  of  structured  grids  is  in  the  dehnition  of  complex 
geometry.  OVERFLOW  2.1  overcomes  this  problem  through  the  use  of  the  Chimera, 
or  overset  grid,  approach.  In  this  approach,  the  grid  system  is  composed  of  overlapping 
structured  grids.  Various  components  of  the  geometry  can  then  be  mapped  to  simple 
structured  grids  and  each  of  the  components  then  compiled  into  one  larger  grid  system. 
Moving  bodies  and  surfaces  can  also  be  easily  simulated  using  overset  grids.  Here,  a 
background  grid  or  set  of  grids  remain  stationary  while  the  grid(s)  associated  with 
the  moving  component  (s)  are  allowed  to  translate  and  rotate.  The  use  of  overset  grids 
in  this  case  precludes  the  need  to  generate  a  new  grid  for  each  time  step  during  the 
motion. 

In  an  overset  grid  system,  information  is  passed  between  grid  blocks  through 
interpolation  of  the  flow  variables.  The  overlapping  points  between  which  information 
is  passed  are  known  as  fringe  points  and  the  interpolation  between  the  fringe  points 
must  be  appropriately  weighted.  In  overset  grids  it  is  also  possible  to  dehne  hole 
points,  wherein  grid  points  from  one  grid  might  not  be  used  in  the  solution  (e.g.  the 
points  of  a  background  grid  which  reside  in  the  volume  of  a  solid  body).  Through 
the  Domain  Connectivity  Function  (DCF),  OVERFLOW  2.1  processes  the  grid  sys¬ 
tem  and  identihes  the  appropriate  hole  cutting,  dehnition  of  fringe  points,  and  grid 
weighting. 
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2-4  Signal  Processing  and  Spectral  Analysis 

In  order  to  examine  the  acoustic  modes  present  in  the  cavity  flow,  it  is  necessary 
to  decompose  the  pressure  time  history  signal  into  its  constituent  frequencies.  This 
decomposition  is  accomplished  through  the  use  of  the  Fourier  transform.  The  classic 
Fourier  transform  is  a  mathematical  operation  which  decomposes  a  signal  by  a  basis 
consisting  of  sines  and  cosines  [19].  For  an  infinite,  continuous  signal  /(t),  the  Fourier 
transform  F{uj)  can  be  calculated  as 


Fiu)  = 


fit) 


(14) 


where  u  represents  frequency,  t  denotes  time,  and  j  is  defined  as  27rh 

By  necessity,  observed  signals  are  finite  and  discrete.  These  properties  require 
the  use  of  the  discrete  Fourier  transform  (DFT)  to  process  these  signals.  The  pro¬ 
cessing  of  finite-duration  signals  leads  to  interesting  behavior  in  the  spectral  analysis. 
Signals  with  frequencies  other  than  those  of  the  basis  set  are  not  periodic  in  the  ob¬ 
servation  window,  which  leads  to  discontinuities  at  the  boundaries  of  the  observation. 
The  result  is  that  these  frequencies  lead  to  non-zero  projections  across  the  entire  basis 
set  [19].  This  phenomenon  is  known  as  spectral  leakage. 

Spectral  leakage  can  be  reduced  through  the  use  of  window  functions.  Win¬ 
dows  are  weighting  functions  applied  to  the  data  set  in  order  to  reduce  the  order  of 
the  discontinuity  at  the  boundary.  The  Hamming  window  is  a  raised-cosine  window 
designed  to  minimize  side  lobe  levels.  This  window  is  defined  as 


27r 

uin)  =  0.54  —  OAGcos  (  —n  )  , 


(15) 


where  N  is  the  length  of  the  discrete  signal  and  n  =  0, 1,  2, . . . ,  iV  —  1.  The  Hamming 
window  is  shown  in  the  time  domain  in  Fig.  5  and  in  the  frequency  domain  in  Fig.  6. 
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Hamming  Window  Frequency  Response  of  Hamming  Window 


Samples 


Figure  5:  The  Hamming  Window  in  Figure  6:  Frequency  Response  of  the 
the  Time  Domain  Hamming  Window 

The  use  of  window  functions  leads  to  a  loss  of  information  near  the  signal  bound¬ 
aries  due  to  the  window  having  small  values  in  these  regions.  In  order  to  avoid  this 
loss  of  data,  overlapping  partitions  can  be  used.  For  a  fixed  signal  length,  50%  over¬ 
lap  in  the  partitions  is  typically  used  in  order  to  hnd  a  balance  between  minimizing 
the  variation  between  the  partitions  and  increased  computational  workload  [20].  The 
Fourier  transform  of  each  windowed  segment  can  then  be  computed  and  the  results 
averaged  to  arrive  at  a  hnal  result.  This  process  of  partitioning  a  signal  into  overlap¬ 
ping  segments,  windowing  each  segment,  computing  the  Fourier  transform  on  each 
segment,  and  then  averaging  the  values  has  come  to  be  known  as  Welch’s  method  [20]. 
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III.  Methodology 

A  CFD  study  using  the  OVERFLOW  2.1  solver  is  conducted  in  order  to  investigate 
the  relationship  between  acoustic  modes  in  a  cavity  flow  and  the  aerodynamic  forces 
and  moments  imposed  on  a  store  released  from  that  cavity.  The  reference  case  chosen 
to  investigate  is  a  cavity  with  L/D  =  4.5,  M  =  0.95,  and  Re  =  2.5  x  This 

cavity  geometry  is  selected  because  it  most  closely  represents  the  typical  geometry  of 
an  aircraft  weapon  bay  of  the  three  geometries  investigated  in  the  WIGS  wind-tunnel 
program.  The  freestream  conditions  correspond  to  Case  3  in  Table  1  and  represent  a 
worst-case  release  point. 

Fonr  cases  are  compared:  the  empty  cavity,  the  cavity  with  a  store  in  the  car¬ 
riage  position,  and  two  configurations  of  the  cavity  with  a  store  in  the  shear  layer. 
These  fonr  cases  are  nominally  depicted  in  Fig.  Ill,  where  case  (a)  corresponds  to  the 
empty  cavity,  (b)  represents  the  cavity  with  the  store  in  carriage  position,  (c)  depicts 
the  cavity  with  the  store  in  the  shear  layer  directly  below  the  carriage  position,  and 
(d)  depicts  the  cavity  with  the  store  in  the  shear  layer  moved  forward  towards  the 
leading  edge  of  the  cavity.  For  convenience,  these  four  situations  will  be  referred  to 
as  the  empty  cavity,  carriage  configuration,  shear  configuration,  and  forward  shear 
configuration,  respectively.  The  inclnsion  of  the  two  shear  configurations  is  to  inves¬ 
tigate  the  finding  from  Heller,  et  ah  [5]  that  a  store  located  in  the  shear  layer  can 
act  as  a  snppression  device  if  it  is  snfficiently  close  to  the  leading  edge  of  the  cavity. 
This  finding  is  of  interest  for  the  current  study  because  if  the  positioning  of  the  store 
as  it  passes  through  the  shear  layer  affects  the  cavity  dynamics,  this  can  affect  the 
trajectory  of  the  released  store. 
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(a)  Empty  cavity 


(b)  Carriage  configuration 


Figure  7:  Depiction  of  four  configurations  investigated 

The  time  pressure  histories  at  six  locations  within  the  cavity  are  recorded  and 
a  spectral  analysis  conducted  on  these  time  histories.  Additionally,  for  the  solutions 
with  a  store  present,  time  pressure  histories  at  nine  locations  on  the  store  and  store 
body  force  and  moment  coefficients  are  calculated  and  a  spectral  analysis  similarly 
conducted  on  these  data.  The  longitudinal  vibrational  modes  of  the  cavity  are  com¬ 
pared  to  the  results  from  the  WIGS  wind  tunnel  test.  The  CFD  derived  store  pressure, 
force,  and  moment  data  are  compared  to  the  predicted  cavity  vibrational  modes  from 
the  CFD  solutions. 


3.1  Grid  Topology  and  Solver  Settings 

The  cavity  itself  is  18  in  long  by  4  in  wide  by  4  in  deep,  matching  the  dimensions 
of  the  WIGS  cavity.  The  computational  domain  extends  approximately  45  in  ahead 
of  the  cavity  and  102  in  behind.  This  domain  is  constructed  of  4  separate  grids  and 
contains  5.2  million  points.  The  cavity  walls  are  constructed  out  of  one  grid  while  the 
interior  of  the  cavity  is  built  out  of  a  second  grid.  This  grid  extends  downward  out 
of  the  cavity  far  enough  to  resolve  the  shear  layer.  The  cavity  grid  is  184  x  45  x  90 
points,  for  a  total  of  745,200  points.  The  plate  which  surrounds  the  cavity,  and  the 
volume  below  this  plate,  are  constructed  from  a  third  grid.  This  grid  wraps  around  the 
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cavity,  leaving  a  hole  in  the  volume  below  the  cavity.  This  empty  volume  necessitates 
the  need  for  a  fourth  grid  to  resolve  the  volume  not  captured  by  the  other  grids.  All 
bay  grids  were  generated  with  a  wall  spacing  of  2  x  10“^  inches,  resulting  in  a  on 
the  order  of  1.  The  cavity  and  store  grids  were  provided  by  Shawn  Westmoreland 
of  Digital  Fusion  Solutions,  Inc.,  and  have  been  demonstrated  to  provide  accurate 
results  in  previous  studies  [3].  Figure  8  shows  the  empty  cavity  grid  system.  The 
cavity  is  represented  by  the  green  area  (1)  and  consists  of  two  grids.  The  blue  grid 

(2)  shows  the  grid  which  defines  the  plate  which  surrounds  the  cavity.  The  red  grid 

(3)  is  the  grid  which  resolves  the  hole  below  the  cavity  which  the  blue  grid  does  not 
dehne.  Figure  9  shows  a  top-down  view  of  the  cavity  grids.  The  green  grid  (1)  dehnes 
the  cavity  walls  and  a  short  lip  around  the  cavity  which  mates  with  the  surrounding 
flat  plate.  The  cavity  itself  is  dehned  by  the  black  grid  (2). 


1 


Figure  8:  Empty  cavity  grid  system 
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Figure  9:  Overhead  view  of  cavity  region 

Six  locations  on  the  cavity  ceiling,  side  wall,  and  aft  wall  are  designated  with 
the  201  boundary  condition  in  OVERFLOW  2.1,  which  allows  for  the  time  history 
of  the  conserved  variables  at  those  locations  to  be  recorded.  Four  of  these  positions 
correspond  to  the  K12,  K16,  K17,  and  K18  pressure  transducers  in  the  WIGS  wind 
tunnel  test.  All  of  these  positions  are  along  the  cavity  centerline.  The  K12  and 
K16  transducers  are  located  on  the  cavity  ceiling  and  the  K17  and  K18  transducers 
are  located  on  the  aft  wall  of  the  cavity.  The  K12  transducer  corresponds  roughly 
to  the  center  of  the  cavity  and  the  K16  transducer  is  nearly  at  the  aft  wall.  The 
K17  transducer  is  roughly  located  at  the  mid-height  of  the  aft  wall  and  the  K18 
transducer  is  approximately  three-quarters  height  off  the  ceiling.  Two  additional 
positions,  LOl  and  L02,  are  located  2  in  down  from  the  cavity  ceiling  on  the  side  wall 
at  longitudinal  positions  which  correspond  to  the  K12  and  K16  transducer  locations, 
respectively.  The  pressures  at  these  locations  were  observed  in  order  to  investigate 
any  lateral  modes  which  may  be  present  in  the  cavity.  The  pressure  probe  locations 
are  illustrated  in  Figs.  10  and  11  and  tabulated  in  Table  3.  The  pressure  histories 
from  these  locations  are  used  in  the  spectral  analysis. 
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K12 


K16 


Figure  10:  Location  of  K12,  K16,  K17,  and  K18  transducers  on  the  cavity  floor 
and  aft  wall 


Figure  11:  Location  of  LOl  and  L02  pressure  probes  on  the  cavity  side  wall 


Table  3:  Location  of  Pressure  Probes 


X  (in) 

y  (in) 

z  (in) 

K12 

9.175 

0.0 

4.0 

K16 

17.725 

0.0 

4.0 

K17 

18.0 

0 

1.975 

K18 

18.0 

0 

0.725 

LOl 

9.175 

2.0 

2.0 

L02 

17.725 

2.0 

2.0 

The  store  of  interest  is  a  1/lOth-scale  modihed  Mk-82  class  weapon.  The  overall 
weapon  shape  is  modelled  but  detailed  features  are  not  included  in  the  model.  This 
is  done  in  order  to  simplify  the  work  required  to  grid  the  store.  Since  the  study  seeks 
to  determine  the  general  relationships  between  the  empty  cavity  flow  and  flow  with 
a  store  present,  and  not  a  store-specihc  solution,  this  compromise  is  acceptable.  The 
surface  and  volume  meshes  for  the  store  are  overly  dense  with  the  longitudinal  grid 
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spacing  never  exceeding  0.075  in.  Altogether,  the  store  grid  system  contains  19  grids 
and  consists  of  4.3  million  grid  points.  The  store  grid  system  is  shown  in  Fig.  12. 


Figure  12:  Store  grid  system 

Three  cases  are  run  with  the  store  present:  the  cavity  with  the  store  in  the 
carriage  position  and  two  conhgurations  with  the  cavity  with  the  store  in  the  shear 
layer.  The  exact  positioning  of  the  store  in  these  three  locations  are  presented  in 
Table  4,  where  the  coordinate  defines  the  location  of  the  nose  of  the  store.  These 
three  situations  are  shown  in  Figs.  13  through  15.  The  flow  direction  in  Figs.  13 
through  15  is  from  left  to  right. 

Table  4:  Location  of  Store  in  Carriage,  Shear,  and  Forward  Shear  Configurations 


X  (in) 

y  (in) 

z  (in) 

Carriage  Configuration 

6.0 

0.0 

2.2 

Shear  Configuration 

6.0 

0.0 

0.0 

Forward  Shear  Conhguration 

2.0 

0.0 

0.0 
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Figure  13:  Cavity  with  the  store  in  carriage  conhguration 


Figure  14:  Cavity  with  the  store  in  shear  conhguration 


Figure  15:  Cavity  with  the  store  in  forward  shear  conhguration 
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Since  the  force  and  moment  data  represent  the  forces  and  moments  acting  on 
the  total  store,  pressure  histories  are  recorded  at  nine  locations  on  the  store  in  order  to 
provide  insight  on  the  contributions  to  the  store  forces  and  moments.  Time  pressure 
histories  are  recorded  at  the  top,  side,  and  bottom  of  the  store  at  each  of  three 
longitudinal  positions:  x/L  =  0.25,  x/L  =  0.5,  and  x/L  =  0.75.  These  positions  are 
shown  in  Fig.  16. 


Figure  16:  Store  pressure  probe  locations 

All  solutions  are  computed  in  OVERFLOW  2.1  unsteady  flow  solver.  OVER¬ 
FLOW  2.1  is  run  using  1st  order  time  with  the  5th  order  spatial  upwind  Mapped 
Weighted  Essentially  Non-Oscillating  (WENOM)  inviscid  flux  algorithm.  A  SSOR 
procedure  is  used  to  solve  the  unfactored  solution  matrix.  The  SST-DDES  turbulence 
model  is  used.  Boundary  conditions  are  established  within  the  domain  to  simulate  a 
leading  plate  length  of  15  in,  corresponding  to  the  WIGS  experimental  setup.  A  time 
step  of  dt  =  1.6  X  10“^s  is  used  in  order  to  capture  the  unsteadiness  of  the  flow  and 
Eve  Newton  subiterations  are  used  between  time  steps.  Nichols  [14]  has  demonstrated 
this  time  step  sufficiently  resolves  the  flow  structures  of  interest  for  this  cavity.  The 
choice  of  Newton  subiterations  is  based  upon  the  result  of  sensitivity  studies,  which 
are  discussed  further  in  the  discussion  of  the  WIGS  cavity  validation.  Grid  resolution 
effects  are  also  investigated. 

The  calculations  are  run  10,000  iterations  and  the  hnal  8,192  time  steps  are 
analyzed.  Time-averaged  data  is  collected  beginning  with  the  1,000th  iteration.  This 
allows  the  simulation  sufficient  time  for  the  transients  associated  with  the  initializa¬ 
tion  of  the  problem  to  be  calculated  out  of  the  solution,  resulting  in  a  physically 
realistic  solution.  At  1,000  iterations,  the  solution  is  statistically  stationary,  which  is 
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shown  in  the  results  in  Chapter  IV.  The  computational  spectral  results  are  computed 
using  Welch’s  method  with  50%  overlap  and  a  Hamming  window  applied  to  each 
data  set.  This  technique  results  in  three  data  sets  of  4,096  samples  each.  The  total 
computational  time  analyzed  is  0.13  s  and  the  time  length  for  each  window  is  0.066 
s.  This  time  period  represents  about  30  cycles  of  the  predominant  460-Hz  frequency 
in  each  window. 

3.2  WICS  Cavity  Validation 

A  solution  for  the  empty  cavity  is  calculated  and  validated  against  the  data  from 
the  WICS  wind  tunnel  test.  This  validation  is  accomplished  in  order  to  establish 
the  proper  solver  settings  which  are  used  in  this  study.  Studies  are  conducted  to 
determine  the  solution’s  sensitivity  to  number  of  Newton  subiterations,  grid  resolution 
and  selection  of  time  step. 

3.2. 1  Newton  Subiteration  Study.  A  Newton  subiteration  study  is  conducted 
to  determine  the  proper  number  of  subiterations  necessary  to  converge  the  solution 
at  each  time  step.  This  study  is  conducted  using  the  medium  grid  and  a  time  step  of 
1.6  X  10“^s.  Solutions  are  computed  using  0,  1,  5,  and  10  Newton  subiterations  and 
are  compared. 

3.2.2  Grid  Resolution  Study.  A  true  grid  independence  study  cannot  be 
performed  for  time-accurate  solutions  which  exhibit  turbulent  structures  because  as 
the  grid  is  rehned,  the  solution  will  resolve  smaller  scale  turbulent  structures.  This 
process  will  continue  until  the  grid  is  refined  below  the  Komologrov  scales,  resulting 
in  a  Direct  Numerical  Simulation  (DNS)  [14].  Grid  resolution  effects  are  studied, 
however,  to  ensure  that  the  solution  resolves  at  least  the  turbulent  structures  of 
interest.  This  study  is  accomplished  by  increasing  and  decreasing  the  rehnement  of 
the  cavity  grid  only.  Adjusting  only  the  cavity  grid  allows  for  the  resolution  of  the 
cavity  to  be  adjusted  without  affecting  the  wall  spacing  on  the  cavity.  The  details  of 
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the  three  grids  used  are  presented  in  Table  5.  The  centerline  plane  of  the  three  grids 
are  shown  in  Fig.  17.  The  grid  resolution  study  was  performed  with  dt  =  1.6  x  10“^. 

Table  5:  Parameters  for  Grid  Rehnement  Study  of  WIGS  Gavity 


Total  Points 

Grid  Dimensions 

^^max  (i^) 

^Umax  (i^) 

AZraax  (In) 

Goarse 

223,200 

124  X  30  X  60 

0.15 

0.15 

0.30 

Medium 

745,200 

184  X  45  X  90 

0.10 

0.10 

0.20 

Fine 

2,049,600 

244  X  70  X  120 

0.08 

0.06 

0.15 

Fine  grid 


Figure  17:  Gomparison  of  three  grid  resolutions  for  cavity  grid 

3.2.3  Time  Step  Sensitivity  Study.  A  time  step  study  is  conducted  using 
the  medium  grid.  Three  solutions  are  run  with  time  steps  of  dt  =  8.0  x  10“®,  dt  = 
1.6  X  10“®,  and  dt  =  4.0  x  10“^  s,  which  correspond  to  time  steps  that  were  investigated 
by  Nichols  [14].  The  three  solutions  are  compared  to  the  WIGS  experimental  data. 


3.3  Cavity  with  Store  in  Carriage  Configuration 

With  the  solution  for  the  empty  cavity  validated,  a  solution  for  the  cavity  with 
the  store  in  the  carriage  conhguration  is  computed.  This  solution  is  run  on  the  medium 
grid  with  dt  =  1.6  x  10“^  s  and  hve  Newton  subiterations.  The  location  of  the  store 
in  the  carriage  conhguration  is  described  above. 

3.4  Cavity  with  Store  in  Shear  Configuration 

One  solution  for  the  cavity  with  the  store  in  the  shear  conhguration  is  computed. 
This  solution  is  run  on  the  medium  grid  with  dt  =  1.6  x  10“^  s  and  hve  Newton 
subiterations.  The  location  of  the  store  in  the  shear  conhguration  is  described  above. 

3.5  Cavity  with  Store  in  Forward  Shear  Configuration 

One  solution  for  the  cavity  with  the  store  in  the  forward  shear  conhguration  is 
computed.  This  solution  is  run  on  the  medium  grid  with  dt  =  1.6  x  10“^  s  and  hve 
Newton  subiterations.  The  location  of  the  store  in  the  forward  shear  conhguration  is 
described  above. 
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IV.  Results 


Several  CFD  simulations  were  conducted  in  order  to  investigate  the  cavity  acoustic 
modes  and  their  relation  to  store  and  moment  loading.  First,  the  simulation  setup  was 
validated  by  running  empty  cavity  solutions  and  comparing  the  results  to  the  WIGS 
experimental  data.  Once  a  set  of  inputs  were  determined  which  provide  an  acceptable 
match  between  the  empty  cavity  simulation  and  WIGS  data,  three  solutions  involving 
the  cavity  and  a  store  were  computed.  For  these  three  solutions,  the  pressure  fluctua¬ 
tions  along  the  cavity  ceiling  and  walls  are  compared  to  the  force  and  moment  loading 
and  pressure  fluctuations  on  the  store.  An  overview  of  the  computational  expense  of 
the  main  solutions  are  presented  in  Table  24  in  Appendix  B.  The  results  of  the  empty 
cavity  validation  and  the  three  solutions  with  the  store  are  presented  next. 

4.1  Empty  WICS  Cavity 

The  empty  cavity  solution  was  hrst  validated  against  the  WIGS  experimental 
data.  In  order  to  determine  the  proper  simulation  parameters,  studies  investigating 
the  solutions’  sensitivity  to  number  of  Newton  subiterations  used,  grid  resolution, 
and  selection  of  time  step  were  conducted.  The  empty  cavity  solution  was  computed 
and  compared  to  the  WIGS  experimental  data.  Additionally,  the  sidewall  pressure 
fluctuations  were  analyzed  and  compared  to  those  on  the  cavity  ceiling  and  aft  wall 
in  order  to  investigate  the  existence  of  any  lateral  modes. 

Newton  Subiteration  Study.  In  order  to  ensure  the  solution  is  fully 
converged  at  each  physical  time  step,  a  Newton  subiteration  study  was  conducted. 
Solutions  using  zero,  one,  hve,  and  ten  Newton  subiterations  were  computed  and 
compared.  In  OVERFLOW  2.1,  specifying  zero  subiterations  calculates  the  solution 
once  for  each  time  step  but  allows  for  the  specihcation  of  using  either  a  first-  or 
second-order  time  accurate  solution.  In  this  study  a  hrst-order  time  accurate  solution 
was  used  due  to  instability  with  the  second-order  accurate  solution.  The  solutions 
were  run  at  a  time  step  of  dt  =  1.6  x  10“^  s.  The  computational  expense  for  each  of 
these  solutions  are  presented  in  Table  25  in  Appendix  B. 


30 


The  SPL  spectra  for  the  four  cases  are  shown  in  Figs.  18  through  21.  The  dashed 
lines  indicate  the  frequencies  of  the  hrst  four  modes  as  predicted  by  the  Rossiter 
equation.  The  WIGS  data  is  plotted  in  the  solid  black  line.  Note  that  WIGS  data 
was  not  available  for  the  K17  location.  The  solution  shows  good  prediction  of  the 
modal  frequencies  with  at  least  hve  Newton  subiterations.  The  solutions  between  hve 
and  ten  Newton  subiterations  do  not  show  a  remarkable  change  in  the  prediction  of 
the  frequencies  nor  amplitudes  of  the  peaks,  specifying  From  these  results,  it  was 
decided  to  use  five  subiterations  for  the  current  research. 


SPL  Spectrum  at  K12  Location 


Figure  18:  Newton  subiteration  sensitivity  for  the  K12  transducer 
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SPL  Spectrum  at  K16  Location 
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Figure  19:  Newton  subiteration  sensitivity  for  the  K16  transducer 
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Figure  20:  Newton  subiteration  sensitivity  for  the  K17  transducer 


32 


SPL  Spectrum  at  K18  Location 


Figure  21:  Newton  subiteration  sensitivity  for  the  K18  transducer 

4.1.2  Grid  Resolution  Study.  Grid  resolution  effects  were  investigated. 
Three  solutions  with  varying  levels  of  rehnement  in  the  cavity  grid  were  computed. 
The  computational  expense  for  each  of  these  solutions  are  presented  in  Table  26  in 
Appendix  B. 

The  sound  pressure  levels  at  the  four  transducer  locations  for  each  of  the  three 
grids  are  presented  in  Figs.  22  through  25.  In  general,  all  three  grids  present  the 
same  solution,  predicting  the  frequency  of  the  same  frequencies  for  the  hrst  four  cavity 
modes,  although  the  amplitudes  do  not  all  agree.  There  are  two  notable  differences 
between  the  grids  in  the  data.  First,  at  the  K12  transducer  location,  the  coarse  grid 
over-predicts  the  overall  sound  levels,  especially  at  the  low  frequencies.  Secondly, 
at  the  K16  transducer  location,  the  hne  grid  under-predicts  the  sound  levels  at  low 
frequencies.  It  is  believed  that  this  is  the  result  of  changing  the  resolution  of  only 
the  cavity  grid  affecting  the  matching  of  cell  sizes  at  the  boundaries  between  grids. 
The  medium  grid  agrees  in  overall  sound  levels  at  all  locations.  From  these  data,  it 
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Sound  Pressure  Level  (dB) 


is  clear  that  all  three  grids  are  rehned  sufficiently  to  resolve  the  turbulent  scales  of 
interest.  Due  to  the  good  agreement  of  sound  levels  across  all  transducer  locations, 
the  medium  grid  provides  the  best  solution  for  the  WIGS  cavity  validation. 


SPL  Spectrum  at  K12  Location 


Figure  22: 


Grid  resolution  sensitivity  for  the  K12  transducer 
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SPL  Spectrum  at  K16  Location 


Figure  23:  Grid  resolution  sensitivity  for  the  K16  transducer 
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Figure  24:  Grid  resolution  sensitivity  for  the  K17  transducer 
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SPL  Spectrum  at  K18  Location 
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Figure  25:  Grid  resolution  sensitivity  for  the  K18  transducer 


4.1.3  Time  Step  Sensitivity  Study.  The  sensitivity  of  the  solution  to  selec¬ 
tion  of  time  step  was  investigated.  Solutions  using  time  steps  of  dt  =  8.0  x  10“®, 
dt  =  1.6  X  10“^,  and  dt  =  4.0  x  10“®  s  were  computed  and  compared.  These  solutions 
were  all  computed  using  hve  Newton  subiterations.  The  computational  expense  for 
each  of  these  solutions  are  presented  in  Table  27  in  Appendix  B. 


The  SPL  spectra  for  the  three  cases  are  shown  in  Figs.  26  through  29.  At  the 
larger  time  step,  the  predicted  frequencies  shift  to  the  left  and  the  solution  predicts 
frequencies  that  are  lower  than  the  experimental  data.  As  the  time  step  is  decreased, 
the  simulation  is  better  able  to  capture  the  periodic  motion  and  the  frequency  is 
better  predicted.  It  is  possible  that  this  solution  could  be  made  more  accurate  if 
more  subiterations  were  used,  although  this  was  not  examined.  The  middle  time 
step  provides  a  good  prediction  of  the  hrst  two  modes  but  slightly  under-predicts 
the  frequency  of  the  higher  modes.  The  smaller  time  step  fully  captures  the  periodic 
motion  and  accurately  predicts  the  frequency  of  all  four  modes,  however  this  solution 


Coarse  grid 
Medium  grid 
Fine  grid 
WICS 
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shows  a  smoothing  of  the  frequency  spectrum.  This  smoothing  is  a  result  of  a  lack 
of  frequency  data  due  to  using  small  data  sets.  Since  each  of  these  three  solutions 
were  run  for  the  same  number  of  time  steps,  the  simulation  physical  time  is  decreased 
as  the  time  step  is  decreased.  The  simulation  with  the  largest  time  step  spans  0.328 
s,  which  includes  151  cycles  at  the  predominant  frequency  of  460  Hz,  whereas  the 
simulation  with  the  smallest  time  step  spans  only  0.066  s,  representing  only  30  cycles 
at  460  Hz.  In  order  to  provide  a  more  robust  solution  with  better-dehned  peaks,  the 
solution  at  the  smallest  time  step  needs  to  be  run  for  a  much  longer  physical  time 
in  order  to  capture  more  cycles  of  the  periodic  motion.  Based  on  a  rule  of  thumb  of 
having  at  least  60  cycles,  the  smallest  time  step  would  need  to  be  run  for  twice  as 
long  as  was  run  here.  It  is  possible  that  some  economy  could  be  gained  from  using  the 
small  time  step,  however,  if  it  results  in  being  able  to  compute  an  accurate  solution 
using  fewer  subiterations,  but  this  was  not  investigated.  Although  the  medium  time 
step  slightly  under-predicts  the  third  and  fourth  modes,  it  provides  an  acceptable 
balance  between  accuracy  and  computational  cost,  and  the  time  step  dt  =  1.6  x  10“® 
s  was  chosen  for  all  of  the  following  solutions. 
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Sound  Pressure  Level  (dB)  Sound  Pressure  Level  (dB) 


SPL  Spectrum  at  K12  Location 


Figure  26:  Time  step  sensitivity  for  the  K12  transducer 


SPL  Spectrum  at  K16  Location 


Figure  27:  Time  step  sensitivity  for  the  K16  transducer 
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Sound  Pressure  Level  (dB)  Sound  Pressure  Level  (dB) 


SPL  Spectrum  at  K17  Location 


Figure  28:  Time  step  sensitivity  for  the  K17  transducer 


SPL  Spectrum  at  K18  Location 


Figure  29:  Time  step  sensitivity  for  the  K18  transducer 
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4-1 -4  Empty  Cavity  Validation.  From  the  sensitivity  studies,  a  set  of  pa¬ 
rameters  were  selected  which  provided  good  agreement  with  the  WIGS  data:  the 
medium  grid,  time  step  of  dt  =  1.6  x  10“^  s,  and  hve  Newton  subiterations.  The 
empty  bay  solution  was  run  with  these  settings  and  the  results  directly  compared  to 
the  WIGS  data.  The  pressure  data  for  the  K18  transducer  is  plotted  against  time  in 
Fig.  30.  All  the  pressure  probe  locations  exhibit  similar  pressure  histories.  The  as¬ 
sumption  set  forth  by  Nichols  [14]  that  demonstrating  that  a  time-accurate  solution  is 
statistically  stationary  is  a  sufficient  condition  for  demonstrating  convergence  is  also 
assumed  here.  It  can  be  shown  that  the  solution  is  statistically  stationary  by  observ¬ 
ing  that  the  statistics  for  each  data  partition  do  not  change.  The  mean  pressure  and 
the  standard  deviation  for  each  data  set  for  each  transducer  location  are  tabulated 
in  Table  6.  The  percentages  below  each  entry  indicates  the  percent  variation  of  that 
parameter  from  the  mean  value  between  data  sets.  None  of  the  mean  pressures  for  a 
data  partition  at  a  given  transducer  location  vary  by  more  than  1.5%  from  the  mean 
for  all  of  the  data.  The  standard  deviations  are  within  14%  for  all  of  the  data  sets 
but  are  within  10%  for  all  but  two  of  the  data  sets.  These  variations  are  sufficiently 
small  such  that  the  statistics  may  be  considered  as  invariant  with  time. 
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Pressure  at  K18  Transducer  Location 


0  0.02  0.04  0.06  0.08  0.1  0.12  0.14  0.16 

Time  (s) 

Figure  30:  Pressure  time  history  for  K18  transducer  location 

The  time  averaged  velocity  vectors  along  the  cavity  centerline  are  shown  in  Fig. 
31.  They  are  colored  by  Mach  number.  The  mean  flow  consists  of  the  shear  layer, 
which  begins  as  the  boundary  layer  at  the  upstream  edge  of  the  cavity  and  spreads 
as  it  travels  downstream.  Upon  approaching  the  aft  wall  of  the  cavity,  part  of  the 
shear  layer  flows  upward  and  out  of  the  cavity,  continuing  downstream  while  the  rest 
stagnates  on  the  aft  wall  and  then  circulates  inside  the  cavity. 
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Table  6:  Mean  pressure  and  standard  deviation  for  each  data  partition  at  each 

transducer  location 


Transducer 

K12 

Mean  pressure  (Pa) 

Standard  deviation  (Pa) 

31,787  31,874  31,875 

-0.09%  0.18%  -0.09% 

1,959  1,596  1,991 

5.99%  -13.68%  7.70% 

Transducer 

K16 

Mean  pressure  (Pa) 

Standard  deviation  (Pa) 

39,174  39,036  39,104 

0.18%  -0.18%  -0.00% 

4,240  3,716  4,251 

4.20%  -8.67%  4.48% 

Transducer 

K17 

Mean  pressure  (Pa) 

Standard  deviation  (Pa) 

36,674  36,591  37,395 

-0.58%  -0.80%  1.38% 

4,129  3,495  4,362 

3.35%  -12.52%  9.18% 

Transducer 

K18 

Mean  pressure  (Pa) 

Standard  deviation  (Pa) 

40,723  40,354  40,647 

0.37%  -0.57%  0.18% 

6,463  5,882  6,092 

5.16%  -4.29%  -0.87% 

The  SPL  spectra  for  the  four  transducer  locations  are  shown  in  Figs.  32  through 
35.  For  the  K12,  K16,  and  K18  locations,  the  WIGS  experimental  data  is  plotted 
alongside  the  CFD  data.  No  WIGS  data  for  the  K17  transducer  location  was  available. 
Again,  the  dashed  lines  indicate  the  frequencies  of  the  hrst  four  modes  as  predicted 
by  the  Rossiter  equation.  All  of  the  results  show  very  good  matching  in  predicting  the 
frequency  of  the  hrst  two  modes,  although  the  amplitude  of  the  second  mode  at  460 
Hz  is  underpredicted  by  approximately  5  dB  at  every  transducer  location.  Note  that 
there  is  no  presence  of  the  hrst  mode  at  the  K12  location  in  either  the  GFD  solution  or 
the  WIGS  experimental  data.  The  third  and  fourth  modes  are  present  but  are  slightly 
misplaced.  The  GFD  solutions  predict  the  third  mode  at  702  Hz,  which  is  30  Hz  lower 
than  the  WIGS  data.  Additionally,  the  GFD  solution  overpredicts  the  amplitude  of 
this  peak  by  approximately  5  dB.  At  all  four  transducer  locations  there  is  a  peak  that 
corresponds  to  the  fourth  mode  at  1025  Hz,  although  in  all  four  locations,  this  peak 
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Figure  31:  Time  averaged  velocity  vectors  for  the  empty  cavity  (colored  by  Mach 
number) 

resides  near  several  other  peaks  of  the  same  or  greater  magnitude,  making  it  unclear 
if  it  truly  represents  the  fourth  mode.  The  K18  location  shows  a  peak  at  976  Hz 
which,  considering  the  downward  shift  of  the  third  mode,  could  possibly  represent  the 
fourth  mode.  The  other  three  locations  have  a  corresponding  peak  at  this  frequency, 
although  it  is  not  as  dominant  in  the  K12,  K16,  nor  K17  locations  as  it  is  in  the  K18 
location.  The  overall  broadband  levels  of  noise  correspond  well  with  the  WIGS  data. 
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SPL  Spectrum  at  K12  Location 


Figure  32:  SPL  spectrum  at  the  K12  transducer  location  for  the  empty  cavity 


SPL  Spectrum  at  K16  Location 


Figure  33:  SPL  spectrum  at  the  K16  transducer  location  for  the  empty  cavity 
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SPL  Spectrum  at  K17  Location 


Figure  34:  SPL  spectrum  at  the  K17  transducer  location  for  the  empty  cavity 


SPL  Spectrum  at  K18  Location 


Figure  35:  SPL  spectrum  at  the  K18  transducer  location  for  the  empty  cavity 
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Pseudo- Schlieren  images  calculated  from  the  pressure  data  along  the  centerline 
of  the  empty  cavity  are  shown  in  Fig.  36.  The  images  are  taken  starting  at  an 
arbitrary  time  t  with  snapshots  at  4.8  x  10“^  s  intervals,  representing  one-quarter  of  a 
cycle  for  the  predominant  cavity  mode  of  460  Hz.  The  pseudo-Schlieren  images  clearly 
show  the  development  of  the  approaching  boundary  layer  into  turbulent  structures 
which  travel  downstream  and  interact  with  the  aft  wall. 


c)  time  =  t  +  9.6e-4  s  d)  time  =  t  +  1.4e-3  s 

Figure  36:  Pseudo-Schlieren  images  of  empty  cavity 
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The  SPL  spectra  for  the  sidewall  locations  are  plotted  in  Figs.  37  and  38.  The 
vertical  dashed  lines  correspond  to  the  frequencies  of  the  hrst  four  longitudinal  modes 
as  predicted  by  the  Rossiter  equation.  These  are  not  plotted  against  any  WICS 
transducer  data  since  the  experiment  did  not  take  any  sidewall  pressure  data.  At 
the  LOl  location,  the  second,  third,  and  fourth  modes  are  clearly  present,  but  the 
hrst  mode  is  indistinguishable.  The  lack  of  a  hrst  mode  coincides  with  the  lack  of  a 
195  Hz  peak  at  the  K12  transducer  location,  which  resides  at  the  same  longitudinal 
distance  in  the  cavity.  At  the  L02  location,  all  four  longitudinal  modes  are  present. 
Similar  trends  are  seen  on  the  sidewall  locations  as  with  the  four  transducer  locations: 
the  second  and  third  modes  show  peaks  at  roughly  the  same  amplitudes,  the  third 
mode  is  predicted  at  702  Hz,  and  the  location  of  the  fourth  mode  is  ambiguous,  with 
peaks  present  at  976  Hz  and  1025  Hz.  In  general,  the  sidewall  pressure  spectra  at 
each  location  matches  the  frequency  content  of  the  transducer  at  the  corresponding 
longitudinal  location,  although  the  amplitude  of  the  signal  on  the  sidewalls  is  reduced 
from  those  on  the  cavity  ceiling. 


Average  SPL  Spectrum  at  L01  Location 


Figure  37:  SPL  spectrum  at  the  LOl  location  for  the  empty  cavity 
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Average  SPL  Spectrum  at  L02  Location 


Figure  38:  SPL  spectrum  at  the  L02  transducer  location  for  the  empty  cavity 

The  turbulence  in  the  cavity  flow  is  dictated  by  the  shedding  of  vortices  from  the 
shear  layer  and  the  fluid  response  to  the  downstream  stagnation  of  these  structures. 
The  feedback  loop  occurs  when  signals  from  the  aft  wall  stagnation  travel  upstream 
and  back  to  the  leading  edge.  This  process  can  essentially  be  mapped  within  the 
volume  of  the  cavity.  As  such,  varying  longitudinal  positions  within  the  cavity  expe¬ 
rience  different  portions  of  the  cavity  flow  phenomenon.  Since  the  longitudinal  modes 
dominate  the  flow,  the  main  direction  of  travel  for  the  signals  are  in  the  streamwise 
directions,  which  explains  why  the  pressure  probes  at  corresponding  longitudinal  lo¬ 
cations  exhibit  similar  frequency  response.  Since  the  sidewalls  essentially  represent 
planes  that  are  parallel  with  the  dominant  flow  mechanism,  the  pressure  fluctuations 
along  the  sidewalls  are  not  as  extreme  as  for  the  ceiling  and  aft  walls,  which  lie  in 
planes  that  are  normal  to  the  flow.  The  lack  of  any  additional  peaks  in  the  sidewall 
pressure  data  indicate  that  there  are  no  signihcant  transverse  waves  traveling  within 
the  cavity.  As  such,  there  is  no  evidence  of  any  lateral  modes,  and  it  is  conhrmed 
here  that  the  longitudinal  modes  dominate  the  flow. 
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The  CFD  solution  is  validated  against  the  WIGS  experimental  data,  show¬ 
ing  good  prediction  of  the  frequency  and  magnitudes  of  the  first  and  second  modes. 
The  solution  shows  reasonable  agreement  with  the  prediction  of  the  third  and  fourth 
modes,  although  the  amplitude  of  the  third  mode  is  overpredicted.  There  is  no  evi¬ 
dence  of  any  lateral  or  vertical  modes  from  the  sidewall  pressure  fluctuation  data.  It 
has  been  shown  that  the  CFD  solution  with  the  settings  shown  here  produce  reason¬ 
able  agreement  with  the  experimental  data. 

4-2  Cavity  with  Store  in  Carriage  Configuration 

With  the  solver  settings  validated  against  the  WIGS  data,  the  solutions  for  the 
carriage,  shear,  and  forward  shear  conhgurations  were  computed.  The  results  for  the 
carriage  conhguration  are  presented  here. 

The  time  averaged  velocity  vectors  along  the  cavity  centerline  are  shown  in  Fig. 
39.  They  are  colored  by  Mach  number.  The  flow  looks  very  similar  to  the  empty 
cavity  scenario.  As  in  the  empty  cavity  flow,  the  shear  layer  spreads  as  it  travels 
downstream,  with  part  of  the  shear  layer  flowing  upward  and  out  of  the  cavity  at 
the  aft  wall  and  the  rest  stagnating  and  then  circulating  inside  the  cavity.  The  flow 
circulates  around  the  cavity,  providing  a  very  similar  prohle  to  what  was  observed 
in  the  empty  cavity,  although  the  flow  between  the  top  of  the  store  and  the  cavity 
ceiling  is  accelerated  due  to  the  geometry  here  being  restricted. 
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Figure  39:  Time  averaged  velocity  vectors  for  the  carriage  configuration  (colored 

by  Mach  number) 


The  SPL  spectra  for  the  four  transducer  locations  are  plotted  against  the  empty 
cavity  results  in  Figs.  40  through  43.  The  vertical  dashed  lines  correspond  to  the 
frequencies  of  the  hrst  four  modes  as  observed  in  the  WIGS  experimental  data.  At 
all  four  transducer  locations,  the  carriage  conhguration  shows  good  agreement  with 
the  empty  cavity  data  in  frequency  and  amplitude  prediction.  The  one  discrepancy 
is  in  the  third  mode,  where  the  carriage  conhguration  shows  a  peak  at  732  Hz,  which 
corresponds  with  the  location  of  the  third  peak  in  the  WIGS  data  and  not  the  third 
peak  that  was  predicted  in  the  empty  cavity  solution.  This  slight  shift  of  one  peak  is 
not  indicative  of  a  store  effect  on  the  acoustic  response.  It  appears  that  the  presence 
of  the  store  in  the  carriage  position  does  not  have  a  signihcant  effect  on  the  huid 
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response  in  the  cavity.  This  observation  was  expected  since  this  was  one  of  the  key 
hndings  from  the  AFFDL  wind  tnnnel  test  [5]. 


SPL  Spectrum  at  K12  Location 


Fignre  40:  SPL  spectrnm  at  the  K12  transdncer  location  for  the  carriage  conhgu- 
ration 
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SPL  Spectrum  at  K16  Location 


Figure  41:  SPL  spectrum  at  the  K16  transducer  location  for  the  carriage  configu¬ 
ration 


SPL  Spectrum  at  K17  Location 
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Figure  42:  SPL  spectrum  at  the  K17  transducer  location  for  the  carriage  conhgu- 
ration 
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SPL  Spectrum  at  K18  Location 


Figure  43:  SPL  spectrum  at  the  K18  transducer  location  for  the  carriage  configu¬ 
ration 

Pseudo- Schlieren  images  from  the  centerline  of  the  cavity  are  shown  in  Fig.  44. 
The  images  are  taken  starting  at  an  arbitrary  time  t  with  snapshots  at  4.8  x  10“^  s 
intervals,  representing  one-quarter  of  a  cycle  for  the  predominant  cavity  mode  of  460 
Hz.  The  presence  of  the  store  does  not  interfere  with  the  vortex  shedding  process 
and  downstream  stagnation.  The  store  is  located  far  enough  up  in  the  cavity  that 
the  shear  layer  passes  below  it.  A  comparison  of  Figs.  36  and  44  shows  a  very  similar 
overall  shape  to  the  structure  in  the  shear  layer. 
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c)  time  =  t  +  9.6e-4  s  d)  time  =  t  +  1.4e-3  s 

Figure  44:  Pseudo-Schlieren  images  of  carriage  configuration 

The  SPL  spectra  for  the  sidewall  locations  are  plotted  in  Figs.  45  and  46.  Again, 
these  results  are  plotted  against  the  spectra  from  the  empty  cavity.  The  same  trends 
which  were  observed  for  the  transducer  locations  are  evident  for  these  cases  as  well, 
including  the  shift  of  the  third  peak  in  the  solution  for  the  carriage  conhguration  and 
the  decreased  overall  sound  levels.  As  with  the  empty  cavity,  the  only  modes  present 
are  the  longitudinal  modes  which  are  observed  along  the  cavity  centerline  and  which 
are  predicted  by  the  Rossiter  equation. 
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SPL  Spectrum  at  L01  Location 
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Figure  45:  SPL  spectrum  at  the  LOl  location  for  the  carriage  conhguration 
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Figure  46:  SPL  spectrum  at  the  L02  transducer  location  for  the  carriage  configu¬ 
ration 
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The  force  and  moment  histories  on  the  store  are  plotted  against  time  in  Figs. 
48  and  49.  The  coordinate  system  is  dehned  as  illnstrated  in  Fig.  47.  The  system 
is  a  right-handed  system  snch  that  the  y-axis  points  into  the  page.  The  mean  and 
standard  deviation  for  each  of  the  forces  and  moments  are  snmmarized  in  Table  7. 
The  dominant  force  on  the  store  is  the  longitndinal  force  imposed  on  the  store.  There 
is  also  a  significant  force  in  the  vertical  direction.  There  is  essentially  no  lateral  force 
on  the  store.  It  is  of  note  that  the  dominant  moment  is  a  pitch-np  moment.  While 
a  pitch-np  behavior  has  been  observed  and  well-docnmented  for  some  cavity-released 
stores,  this  behavior  is  often  attribnted  to  the  interaction  of  the  store  with  the  shear 
layer.  These  data  snggest,  however,  that  the  store  is  snbjected  to  a  pitch-np  moment 
from  the  time  of  release.  There  is  also  a  flnctnating  roll  moment,  althongh  this  is 
centered  abont  the  x-axis.  There  is  effectively  no  yaw  moment  present  on  the  store. 


Z 


Fignre  47:  Dehnition  of  coordinate  system 
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Figure  48:  Force  time  history  for  the  store  in  the  carriage  conhguration 
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Figure  49:  Moment  time  history  for  the  store  in  the  carriage  conhguration 
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Table  7:  Summary  of  store  forces  and  moments  in  carriage  configuration 


Mean  Value 

Standard  Deviation 

Longitudinal  force 

14.00  N 

9.33  N 

Lateral  force 

0.01  N 

0.01  N 

Vertical  force 

-0.55  N 

2.23  N 

Roll  moment 

-0.04  N-m 

0.12  N-m 

Pitch  moment 

1.16  N-m 

2.98  N-m 

Yaw  moment 

0.00  N-m 

0.01  N-m 

The  time-averaged  pressure  coefficient  along  the  top  and  bottom  of  the  store 
are  presented  in  Fig.  50.  There  is  an  overwhelming  negative  pressure  acting  on  the 
majority  of  the  store.  The  pressure  along  the  top  of  the  store  steadily  increases  from 
xj L  =  0.0  io  x/L  =  0.7.  Over  this  same  span,  the  pressure  along  the  bottom  of  the 
store  oscillates  but  displays  a  slight  overall  decrease  in  pressure.  Aft  oi  x/L  =  0.7, 
interference  with  the  fins  on  the  store  creates  a  sharp  increase  in  pressure,  although 
there  maintains  a  pressure  differential  between  the  top  and  bottom.  From  the  pressure 
distribution,  the  center  of  pressure  can  be  calculated.  Table  8  summarizes  the  average 
pressure  coefficient  differential  (taken  as  the  pressure  coefficient  on  the  top  of  the  store 
minus  the  pressure  coefficient  on  the  bottom  of  the  store),  the  center  of  pressure,  and 
the  store’s  center  of  mass.  From  the  dehnition  of  the  pressure  coefficient  differential 
used  here,  a  positive  value  indicates  a  net  force  downward  on  the  store.  The  pitch-up 
behavior  can  thus  be  explained:  the  net  pressure  differential  creates  a  downward  force 
which  acts  at  a  point  on  the  store  which  is  aft  of  the  center  of  mass,  resulting  in  a 
pitch- up  moment.  Note  that  with  the  center  of  pressure  aft  of  the  center  of  gravity, 
the  store  is  stable. 
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Figure  50:  Pressure  coefficient  on  top  and  bottom  of  store  in  carriage  configuration 


Table  8:  Average  pressure  differential,  center  of  pressure,  and  center  of  mass  for 

the  store  in  carriage  configuration 


Average  Pressure  Differential 
Cp^top  C 

Center  of  Pressure 
x/L 

Center  of  Mass 
x/L 

0.041 

0.569 

0.459 

The  longitudinal  and  vertical  force  spectra  are  plotted  alongside  the  K16  trans¬ 
ducer  spectrum  in  Figs.  51  and  52.  The  K16  transducer  was  selected  because  it 
showed  the  best  prediction  of  the  hrst  four  acoustic  modes.  Note  that  the  longitudi¬ 
nal  and  vertical  forces  are  plotted  on  different  scales  and  differ  by  about  two  orders  of 
magnitude.  The  force  spectra  are  plotted  on  a  logarithmic  scale  in  order  to  match  the 
scaling  with  the  cavity  ceiling  SPL.  The  force  spectra  show  a  very  strong  correlation 
to  the  pressure  fluctuations  on  the  cavity  ceiling.  The  hrst  four  modes  as  observed 
on  the  cavity  ceiling  are  also  apparent  in  the  store  forces.  The  relative  amplitudes  of 
the  peaks  from  the  broadband  noise  levels  also  correlate.  The  pitch  and  roll  moment 
spectra  are  also  plotted  alongside  the  K16  transducer  spectrum  in  Figs.  53  and  54. 
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As  with  the  force  spectra,  the  moment  spectra  are  plotted  on  a  logarithmic  scale.  As 
with  the  forces,  note  that  the  pitch  and  roll  moments  are  plotted  on  different  scales 
and  differ  by  about  three  orders  of  magnitude.  The  pitch  moment  spectrum  shows 
a  very  strong  correlation  to  the  K16  SPL  spectrum.  The  modal  frequencies  and  the 
peak  relative  amplitudes  match  very  well  with  the  K16  SPL.  The  roll  moment  spec¬ 
trum  does  not  show  any  correlation  to  the  K16  spectrum.  None  of  the  cavity  modes 
are  present  in  the  roll  moment  data.  There  is  one  peak  in  the  roll  moment  spectrum 
which  corresponds  to  the  third  mode,  but  its  amplitude  does  not  distinguish  it  from 
the  broadband  noise. 


Cavity  Ceiling  Pressure  and  Store  Longitudinal  Force  Spectra 


Figure  51:  Store  longitudinal  force  frequency  spectrum  (right  axis)  and  K16  trans¬ 
ducer  SPL  spectrum  (left  axis)  for  carriage  conhguration 
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Cavity  Ceiling  Pressure  and  Store  Vertical  Force  Spectra 


Figure  52:  Store  vertical  force  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  carriage  configuration 


Cavity  Ceiling  Pressure  and  Store  Pitch  Moment  Spectra 


Figure  53:  Store  pitch  moment  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  carriage  conhguration 
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Cavity  Ceiling  Pressure  and  Store  Roll  Moment  Spectra 
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Figure  54:  Store  roll  moment  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  carriage  configuration 


The  SPL  spectra  on  the  surface  of  the  store  are  presented  in  Figs.  55  through 
57.  Each  of  the  hgures  presents  the  SPL  spectra  at  a  given  longitudinal  position  on 
the  store.  The  pressure  spectra  at  each  longitudinal  position  correlate  well  with  each 
other.  At  x/L  =  0.25,  the  second,  third,  and  fourth  modes  are  clearly  present  with  the 
second  mode  being  the  dominant  frequency.  This  longitudinal  position  corresponds 
almost  with  the  longitudinal  position  of  the  K12  transducer,  so  it  is  not  surprising  to 
see  a  lack  of  the  hrst  mode  here.  At  x/L  =  0.50,  the  hrst  three  modes  are  clearly  evi¬ 
dent,  although  here  the  second  and  third  modes  are  roughly  equivalent  in  amplitude. 
At  the  x/L  =  0.75  location,  the  only  clear  mode  present  is  the  third  mode. 
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Average  Pressure  at  X/L  =  0.25 


Figure  55:  SPL  spectra  on  store  body  at  x/L  =  0.25  for  carriage  configuration 


Average  Pressure  at  X/L  =  0.50 


Figure  56:  SPL  spectra  on  store  body  at  x/L  =  0.50  for  carriage  configuration 
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Average  Pressure  at  X/L  =  0.75 


Figure  57:  SPL  spectra  on  store  body  at  x/L  =  0.75  for  carriage  configuration 

The  same  vibrational  modes  which  are  observed  in  the  empty  cavity  are  also 
present  for  the  cavity  with  a  store  in  carriage.  Previous  observations  that  small  (rela¬ 
tive  to  the  cavity)  stores  located  in  the  cavity  do  not  affect  the  acoustic  response  of  the 
cavity  are  conhrmed  here.  A  comparison  of  the  cavity  surface  pressure  fluctuations 
and  the  store  force  and  moment  loadings  shows  correlation  between  the  longitudi¬ 
nal  force,  vertical  force,  and  pitch  moment  with  the  cavity  vibrational  modes.  The 
dominant  moment  on  the  store  in  the  carriage  location  is  a  pitch-up  moment.  The 
pressure  fluctuations  on  the  store  in  the  cavity  vary  by  longitudinal  position  and 
correlate  between  circumferential  position  at  a  given  longitudinal  position. 

4.3  Cavity  with  Store  in  Shear  Configuration 

The  results  for  the  shear  configuration  are  presented  here.  It  is  shown  that 
the  flow  characteristics  in  the  shear  configuration  differ  greatly  from  the  carriage 
conhguration.  It  is  expected  that  these  fundamental  differences  will  influence  the 
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behavior  of  a  cavity-released  store  as  it  transitions  between  the  two  positions  during 
a  separation  event. 

The  time  averaged  velocity  vectors  along  the  cavity  centerline  are  shown  in 
Fig.  58.  They  are  colored  by  Mach  number.  The  presence  of  the  store  prevents  the 
shear  layer  from  from  spreading  as  in  the  empty  cavity  and  the  carriage  conhguration. 
Instead,  the  shear  layer  splits  as  two  thinner  layers  travel  around  the  top  and  bottom 
of  the  store,  reforming  behind  the  store.  Once  it  has  reformed,  however,  it  behaves 
similarly  as  in  the  previous  conhgurations,  with  part  of  the  shear  layer  flowing  upward 
and  out  of  the  cavity  at  the  aft  wall  and  the  rest  stagnating  and  then  circulating  inside 
the  cavity. 
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Figure  58:  Time  averaged  velocity  vectors  for  the  shear  conhguration  (colored  by 
Mach  number) 
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The  SPL  spectra  for  the  four  transducer  locations  are  plotted  against  the  empty 
cavity  results  in  Figs.  59  through  62.  The  vertical  dashed  lines  correspond  to  the 
frequencies  of  the  hrst  four  modes  as  predicted  by  the  Rossiter  equation.  Generally 
speaking,  the  shear  conhguration  shows  good  agreement  with  the  empty  cavity  data 
in  frequency  and  amplitude  prediction  at  all  four  transducer  locations.  As  with  the 
carriage  conhguration,  the  shear  conhguration  shows  a  peak  at  732  Hz  for  the  third 
mode,  which  corresponds  with  the  location  of  the  third  peak  in  the  WIGS  data. 
Also,  the  fourth  mode  at  the  K12  location  is  much  more  clearly  dehned  in  the  shear 
conhguration  than  the  empty  cavity.  Additionally,  the  overall  noise  levels  and  the 
amplitude  of  the  hrst  mode  are  lower  in  the  shear  conhguration  than  the  empty 
cavity.  This  slight  suppression  of  sound  levels  is  possibly  the  result  of  the  store  acting 
as  a  spoiler. 


SPL  Spectrum  at  K12  Location 


Figure  59:  SPL  spectrum  at  the  K12  transducer  location  for  the  shear  conhguration 
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SPL  Spectrum  at  K16  Location 


Figure  60: 


Figure  61: 


SPL  spectrum  at  the  K16  transducer  location  for  the  shear  conhguration 


SPL  Spectrum  at  K17  Location 
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SPL  spectrum  at  the  K17  transducer  location  for  the  shear  conhguration 
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SPL  Spectrum  at  K18  Location 


Figure  62:  SPL  spectrum  at  the  K18  transducer  location  for  the  shear  conhguration 

Pseudo- Schlieren  images  from  the  centerline  of  the  empty  cavity  are  shown  in 
Fig.  63.  The  images  are  taken  starting  at  an  arbitrary  time  t  with  snapshots  at 
4.8  X  10“^  s  intervals,  representing  one-quarter  of  a  cycle  for  the  predominant  cavity 
mode  of  460  Hz.  Here  the  store  clearly  interacts  with  the  shear  layer.  The  shear 
layer  is  forced  to  flow  around  the  store,  although  the  vortical  structures  generally 
remain  intact  and  continue  to  flow  downstream  where  they  eventually  stagnate  at 
the  aft  cavity  wall.  The  vortex  shedding/cavity  response  cycle  is  largely  unimpeded 
despite  the  presence  of  the  store  in  the  shear  layer.  This  observation  supports  the 
findings  that  the  pressure  spectra  are  relatively  unchanged  from  the  empty  cavity 
conhguration.  There  does  appear  to  be  some  dehection  of  the  shear  layer  in  snapshot 
(d),  wherein  the  store  causes  the  shear  layer  to  how  down  and  out  of  the  cavity  on  an 
intermittent  basis.  This  dehection  of  the  how  would  explain  the  slightly  lower  cavity 
noise  levels  observed  at  the  K18  transducer  at  low  frequencies. 
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c)  time  =  t  +  9.6e-4  s  d)  time  =  t  +  1.4e-3  s 

Figure  63:  Pseudo-Schlieren  images  of  shear  configuration 


The  SPL  spectra  for  the  sidewall  locations  are  plotted  in  Figs.  64  and  65.  Again, 
these  results  are  plotted  against  the  spectra  from  the  empty  cavity.  As  with  the  four 
transducer  locations,  these  spectra  show  good  correlation  with  the  empty  cavity  and 
exhibit  the  same  trends.  Namely,  the  third  mode  for  the  shear  configuration  better 
corresponds  to  the  observed  frequency  from  the  WIGS  experimental  data  than  the 
empty  cavity,  the  fourth  mode  at  the  LOl  location  is  more  clearly  dehned  in  the  shear 
conhguration,  and  the  amplitude  of  the  hrst  peak  at  the  L02  location  is  predicted  to 
be  lower  in  the  shear  conhguration  than  the  empty  cavity. 
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Figure  64:  SPL  spectrum  at  the  LOl  location  for  the  shear  configuration 


SPL  Spectrum  at  L01  Location 
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Figure  65:  SPL  spectrum  at  the  L02  transducer  location  for  the  shear  conhguration 
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The  force  and  moment  histories  on  the  store  are  plotted  against  time  in  Figs. 
66  and  67.  These  time  histories  exhibit  the  same  general  characteristics  as  those 
for  the  carriage  conhgnration.  The  mean  and  standard  deviation  for  each  of  the 
forces  and  moments  are  snmmarized  in  Table  9.  The  dominant  force  on  the  store  is 
the  longitndinal  force  imposed  on  the  store.  There  is  also  a  signihcant  force  in  the 
vertical  direction.  There  is  essentially  no  lateral  force  on  the  store.  The  dominant 
moment  is  a  pitch  np  moment  and  there  is  also  a  flnctnating  roll  moment,  althongh 
this  is  centered  abont  the  x-axis.  There  is  effectively  no  yaw  moment  present  on  the 
store. 
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Force  time  history  for  the  store  in  the  shear  conhgnration 


71 


Moments 


0.08  0.1 
Time  (s) 


0.16 


Figure  67:  Moment  time  history  for  the  store  in  the  shear  conhguration 


Table  9:  Summary  of  store  forces  and  moments  in  shear  conhguration 


Mean  Value 

Standard  Deviation 

Longitudinal  force 

11.19  N 

7.67  N 

Lateral  force 

0.00  N 

0.01  N 

Vertical  force 

-0.37  N 

2.52  N 

Roll  moment 

-0.02  N-m 

0.10  N-m 

Pitch  moment 

0.73  N-m 

2.87  N-m 

Yaw  moment 

0.00  N-m 

0.01  N-m 

The  time-averaged  pressure  coefficient  along  the  top  and  bottom  of  the  store  are 
presented  in  Fig.  68.  The  pressures  at  the  top  and  bottom  of  the  store  are  almost  equal 
along  the  entire  length  of  the  store.  This  distribution  is  a  consequence  of  the  symmetry 
of  the  store  and  both  surfaces  essentially  seeing  the  same  how.  There  is  a  high  pressure 
at  the  nose  of  the  store  from  the  approaching  shear  layer,  but  the  pressures  over  most 
of  the  store  are  almost  equal  to  the  freestream.  This  is  markedly  diherent  from  the 
pressure  distribution  seen  on  the  store  in  the  carriage  position.  Table  10  summarizes 
the  average  pressure  coefficient  diherential,  the  center  of  pressure,  and  the  store’s 
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center  of  mass.  As  with  the  store  in  the  carriage  conhguration,  a  pitch-up  moment 
is  observed,  although  the  mechanism  here  is  different.  There  is  a  net  upward  force 
acting  at  a  point  ahead  of  the  center  of  mass  that  creates  the  pitch-up  moment.  Note 
that  with  the  center  of  pressure  ahead  of  the  center  of  gravity,  the  store  is  unstable 
here. 


Coefficient  of  Pressure 


Figure  68:  Pressure  coefficient  on  top  and  bottom  of  store  in  shear  conhguration 


Table  10:  Average  pressure  differential,  center  of  pressure,  and  center  of  mass  for 
the  store  in  shear  conhguration 


Average  Pressure  Differential 

Cp^top  C p^}jottom 

Center  of  Pressure 
x/L 

Center  of  Mass 
x/L 

-0.047 

-0.078 

0.459 

The  mean  and  standard  deviation  of  the  longitudinal  force  and  pitch  moment 
between  the  carriage  and  shear  conhgurations  are  compared  in  Table  11.  The  forces 
and  moments  on  the  store  are  greater  in  the  carriage  position  than  in  the  shear  layer. 
The  carriage  conhguration  data  is  also  more  variant,  indicating  a  wider  range  of 
huctuations.  This  observation  initially  seems  counter-intuitive  since  greater  pitch-up 
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behavior  has  typically  been  observed  to  occur  as  a  released  store  approaches  the  shear 


layer.  However,  this  can  be  explained  when  the  store’s  stability  is  taken  into  consid¬ 
eration.  As  the  store  travels  downward  in  the  cavity,  it  transitions  from  a  stable  to  an 
unstable  conhguration  as  the  center  of  pressure  moves  forward.  Although  the  store 
experiences  higher  moments  early  in  the  trajectory,  its  stable  conhguration  makes 
it  more  resistant  to  pitching.  As  the  store  transitions  to  an  unstable  conhguration, 
however,  the  store  is  less  able  to  self-correct  and  the  pitching  moment  may  be  more 
readily  observed  in  the  trajectory  of  the  store. 


Table  11:  Comparison  of  store  forces  and  moments  between  carriage  and  shear 

conhgurations 


Carriage  conhguration 

Shear  conhguration 

Mean 

Standard  Deviation 

Mean 

Standard  Deviation 

Longitudinal 

force 

14.00  N 

9.33  N 

11.19  N 

7.67  N 

Vertical 

force 

-0.55  N 

2.23  N 

-0.37  N 

2.52  N 

Pitch 

moment 

1.16  N-m 

2.98  N-m 

0.73  N-m 

2.87  N-m 

The  longitudinal  and  vertical  force  spectra  are  plotted  alongside  the  K16  trans¬ 
ducer  spectrum  in  Figs.  69  and  70.  The  force  spectra  are  plotted  on  a  logarithmic 
scale  in  order  to  match  the  scaling  with  the  cavity  ceiling  SPL.  The  force  spectra 
show  a  strong  correlation  to  the  pressure  huctuations  on  the  cavity  ceiling  and  simi¬ 
lar  trends  that  were  observed  in  the  carriage  conhguration  are  present  here,  although 
the  correlation  with  the  longitudinal  force  spectrum  is  not  as  strong  here  as  in  the 
carriage  conhguration.  The  hrst  four  modes  as  observed  on  the  cavity  ceiling  are  also 
apparent  in  the  store  forces.  The  relative  amplitudes  of  the  peaks  from  the  broad¬ 
band  noise  levels  also  correlate.  The  pitch  and  roll  moment  spectra  are  also  plotted 
alongside  the  K16  transducer  spectrum  in  Figs.  71  and  72.  As  with  the  force  spectra, 
the  moment  spectra  are  plotted  on  a  logarithmic  scale.  Again,  note  that  the  pitch 
and  roll  moments  are  plotted  on  diherent  scales  and  diher  by  about  three  orders  of 
magnitude.  The  pitch  moment  spectrum  shows  a  very  strong  correlation  to  the  K16 
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SPL  spectrum.  The  modal  frequencies  and  the  peak  relative  amplitudes  match  very 
well  with  the  K16  SPL.  The  roll  moment  spectrum  does  not  show  any  correlation  to 
the  K16  spectrum.  None  of  the  cavity  modes  are  present  in  the  roll  moment  data. 


Cavity  Ceiling  Pressure  and  Store  Longitudinal  Force  Spectra 


Figure  69:  Store  longitudinal  force  frequency  spectrum  (right  axis)  and  K16  trans¬ 
ducer  SPL  spectrum  (left  axis)  for  shear  configuration 
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Cavity  Ceiling  Pressure  and  Store  Vertical  Force  Spectra 
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Figure  70:  Store  vertical  force  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  shear  conhguration 


Cavity  Ceiling  Pressure  and  Store  Pitch  Moment  Spectra 


Figure  71:  Store  pitch  moment  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  shear  conhguration 
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Cavity  Ceiling  Pressure  and  Store  Roll  Moment  Spectra 


Figure  72:  Store  roll  moment  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  shear  conhguration 

The  SPL  spectra  on  the  surface  of  the  store  are  presented  in  Figs.  73  through 
75.  Each  of  the  hgures  presents  the  SPL  spectra  at  a  given  longitudinal  position  on 
the  store.  The  pressure  spectra  at  each  longitudinal  position  correlate  well  with  each 
other.  At  x/L  =  0.25,  the  second,  third,  and  fourth  modes  are  clearly  present  with 
the  second  mode  being  the  dominant  frequency.  The  hrst  mode  is  present,  although 
it  is  located  at  122  Hz,  which  is  lower  than  the  hrst  mode  observed  on  the  cavity 
ceiling.  At  x/L  =  0.50,  the  second  and  third  modes  are  clearly  evident,  although 
here  the  second  and  third  modes  are  roughly  equivalent  in  amplitude.  At  the  x/L 
=  0.75  location,  the  hrst,  third,  and  fourth  modes  are  present.  As  with  the  x/L  = 
0.25  location,  the  hrst  mode  is  predicted  at  a  frequency  lower  than  observed  on  the 
cavity  ceiling.  These  are  essentially  the  same  trends  that  were  observed  in  the  carriage 
conhguration.  While  the  characteristics  of  the  mean  how  are  seen  to  diher  between 
the  two  positions,  the  characteristics  of  the  turbulence  in  these  locations  correspond 
very  well,  as  they  are  driven  by  the  vortex  shedding  and  cavity  huid  response  which 
dehnes  the  how  throughout  the  entire  cavity. 
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Average  Pressure  at  X/L  =  0.25 


Figure  73:  SPL  spectra  on  store  body  at  x/L  =  0.25  for  shear  configuration 


Average  Pressure  at  X/L  =  0.50 


Figure  74:  SPL  spectra  on  store  body  at  x/L  =  0.50  for  shear  configuration 
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Average  Pressure  at  X/L  =  0.75 


Figure  75:  SPL  spectra  on  store  body  at  x/L  =  0.75  for  shear  configuration 

The  same  vibrational  modes  which  are  observed  in  the  empty  cavity  are  also 
present  for  the  cavity  with  a  store  in  the  shear  layer.  Although  previous  experiments 
have  shown  that  a  store  located  close  to  the  cavity  leading  edge  can  act  as  a  suppressor, 
it  is  shown  here  that  if  the  weapon  is  sufficiently  far  from  the  cavity  leading  edge,  it 
will  not  have  an  effect  on  the  acoustic  response  of  the  cavity.  A  comparison  of  the 
cavity  surface  pressure  ffuctuations  and  the  store  force  and  moment  loadings  shows 
correlation  between  the  longitudinal  force,  vertical  force,  and  pitch  moment  with  the 
cavity  vibrational  modes.  The  dominant  moment  on  the  store  in  the  shear  location 
is  a  pitch- up  moment.  The  pressure  fluctuations  on  the  store  in  the  cavity  vary 
by  longitudinal  position  and  correlate  between  circumferential  position  at  a  given 
longitudinal  position.  The  cavity  response  and  store  force/moment  loading  for  this 
conhguration  are  consistent  with  what  was  observed  for  the  carriage  conhguration. 
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4-4  Cavity  with  Store  in  Forward  Shear  Configuration 

The  results  for  the  forward  shear  conhguration  are  presented  here.  The  over¬ 
all  characteristics  of  the  flow  are  not  greatly  different  from  the  two  configurations 
presented  earlier,  although  there  are  notable  differences. 

The  time  averaged  velocity  vectors  along  the  cavity  centerline  are  shown  in  Fig. 
76.  They  are  colored  by  Mach  number.  As  with  the  shear  conhguration,  the  presence 
of  the  store  forces  the  shear  layer  to  split  as  the  how  travels  around  the  top  and 
bottom  of  the  store,  reforming  behind  the  store.  Unlike  in  the  shear  conhguration, 
however,  the  store  is  upstream  far  enough  such  that  the  shear  layer  is  still  able  to 
spread  out  before  reaching  the  aft  cavity  wall.  In  this  way,  the  time-averaged  how 
much  resembles  that  of  the  empty  cavity  and  the  carriage  conhguration. 


Figure  76:  Time  averaged  velocity  vectors  for  the  forward  shear  conhguration  (col¬ 
ored  by  Mach  number) 
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The  SPL  spectra  for  the  four  transducer  locations  are  plotted  against  the  empty 
cavity  results  in  Figs.  77  through  80.  The  vertical  dashed  lines  correspond  to  the 
frequencies  of  the  hrst  four  modes  as  predicted  by  the  Rossiter  equation.  Generally 
speaking,  the  forward  shear  conhguration  shows  good  agreement  with  the  empty 
cavity  data  in  frequency  and  amplitude  prediction  at  all  four  transducer  locations. 
As  with  the  carriage  and  shear  configurations,  the  forward  shear  configuration  shows 
a  peak  at  732  Hz  for  the  third  mode,  which  corresponds  with  the  location  of  the  third 
peak  in  the  WIGS  data.  This  conhguration  shows  better  matching  with  the  WIGS 
empty  cavity  data  in  the  amplitude  of  the  second  and  third  modes  at  all  transducer 
locations;  specihcally  the  amplitude  of  the  second  mode  is  predicted  higher  and  the 
amplitude  of  the  third  mode  is  predicted  lower  than  the  empty  cavity  solution.  One 
signihcant  difference  is  that  the  amplitude  of  the  fourth  mode  is  greatly  increased  at 
the  K12,  K16,  and  K17  locations.  Also,  in  the  K12  location,  there  is  now  a  distinct 
peak  at  the  hrst  mode.  The  overall  noise  levels  are  in  agreement  with  the  empty 
cavity  solution  across  all  transducer  locations,  which  indicates  that  the  store  is  not 
acting  as  a  suppressor  in  this  situation. 
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SPL  Spectrum  at  K12  Location 


Figure  77:  SPL  spectrum  at  the  K12  transducer  location  for  the  forward  shear 

conhguration 


SPL  Spectrum  at  K16  Location 


Figure  78:  SPL  spectrum  at  the  K16  transducer  location  for  the  forward  shear 

conhguration 
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SPL  Spectrum  at  K17  Location 
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Figure  79:  SPL  spectrum  at  the  K17  transducer  location  for  the  forward  shear 

conhguration 
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SPL  Spectrum  at  K18  Location 


Figure  80:  SPL  spectrum  at  the  K18  transducer  location  for  the  forward  shear 

conhguration 
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Pseudo-Schlieren  images  from  the  centerline  of  the  empty  cavity  are  shown  in 
Fig.  81.  The  images  are  taken  starting  at  an  arbitrary  time  t  with  snapshots  at 
4.8  X  10“"^  s  intervals,  representing  one-quarter  of  a  cycle  for  the  predominant  cavity 
mode  of  460  Hz.  Here,  as  with  the  shear  conhguration,  the  store  clearly  interacts 
with  the  shear  layer.  The  shear  layer  is  forced  to  flow  around  the  store,  although  the 
vortical  structures  generally  remain  intact  and  continue  to  flow  downstream  where 
they  eventually  stagnate  at  the  aft  cavity  wall.  With  the  store  located  sufficiently 
upstream,  there  is  ample  time  for  the  shear  layer  to  reform  behind  the  store,  and  the 
vortex  shedding/cavity  response  cycle  is  largely  unimpeded  despite  the  presence  of  the 
store  in  the  shear  layer.  This  further  supports  the  hndings  that  the  pressure  spectra 
are  relatively  unchanged  from  the  empty  cavity  conhguration.  It  also  fails  to  support 
earlier  observations  that  a  store  located  in  the  shear  layer  acts  to  suppress  the  vortex 
shedding/acoustic  response  cycle  [5].  The  stores  used  in  the  AFFDL  experiment  were 
cylindrical  with  an  ogive  nosecone  and  were  longer  relative  to  the  overall  cavity  length 
than  the  store  model  used  here.  Furthermore,  the  exact  conhgurations  for  which  a 
suppressive  effect  was  observed,  including  positioning  of  the  store,  are  not  detailed. 
Just  as  different  spoiler  designs  have  been  shown  to  have  different  effectiveness  in 
the  suppression  of  vibrations,  the  exact  store  shape,  size,  and  positioning  would  be 
expected  to  have  differing  effects  on  the  huid  cavity  response  studied  here. 
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c)  time  =  t  +  9.6e-4  s  d)  time  =  t  +  1.4e-3  s 


Figure  81:  Pseudo-Schlieren  images  of  forward  shear  configuration 

The  SPL  spectra  for  the  sidewall  locations  are  plotted  in  Figs.  82  and  83.  Again, 
these  results  are  plotted  against  the  spectra  from  the  empty  cavity.  Similar  trends 
as  were  observed  with  the  four  transducer  locations  are  observed  in  these  spectra. 
Namely,  the  amplitude  of  the  second  mode  is  predicted  higher  than  the  empty  cavity, 
the  amplitude  of  the  third  mode  is  predicted  lower  than  the  empty  cavity,  and  the 
amplitude  of  the  fourth  mode  is  predicted  to  increase  from  the  empty  cavity  solution. 
As  with  the  empty  cavity  solution,  there  appears  to  be  no  discernable  hrst  mode  at 
the  L02  location. 
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Figure  82:  SPL  spectrum  at  the  LOl  location  for  the  forward  shear  conhguration 
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Figure  83:  SPL  spectrum  at  the  L02  transducer  location  for  the  forward  shear 

configuration 


The  force  and  moment  histories  on  the  store  are  plotted  against  time  in  Figs. 
84  and  85.  These  time  histories  exhibit  the  same  general  characteristics  as  those  for 
the  carriage  and  shear  conhgnrations.  The  mean  and  standard  deviation  for  each  of 
the  forces  and  moments  are  snmmarized  in  Table  12.  The  dominant  force  on  the  store 
is  the  longitndinal  force  imposed  on  the  store.  There  is  also  a  signihcant  force  in  the 
vertical  direction.  There  is  essentially  no  lateral  force  on  the  store.  The  dominant 
moment  is  a  pitch  np  moment  and  there  is  also  a  flnctnating  roll  moment,  althongh 
this  is  centered  abont  the  x-axis.  There  is  effectively  no  yaw  moment  present  on  the 
store. 
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Fignre  84:  Force  time  history  for  the  store  in  the  forward  shear  conhgnration 
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Figure  85: 

Table  12: 


Mean  Value 

Standard  Deviation 

Longitudinal  force 

11.64  N 

9.01  N 

Lateral  force 

0.00  N 

0.01  N 

Vertical  force 

-0.38  N 

2.68  N 

Roll  moment 

-0.02  N-m 

0.10  N-m 

Pitch  moment 

1.14  N-m 

3.17  N-m 

Yaw  moment 

0.00  N-m 

0.01  N-m 

Moment  time  history  for  the  store  in  the  forward  shear  conhguration 
Summary  of  store  forces  and  moments  in  forward  shear  configuration 


The  time-averaged  pressure  coefficient  along  the  top  and  bottom  of  the  store 
are  presented  in  Fig.  86.  It  looks  remarkably  similar  to  the  shear  conhguration. 
The  pressures  at  the  top  and  bottom  of  the  store  are  almost  equal  along  the  entire 
length  of  the  store.  The  pressure  coefficient  along  the  bottom  of  the  store  shows  the 
same  huctuations  at  the  same  longitudinal  positions  on  the  store  as  the  solution  from 
the  shear  conhguration.  One  distinguishing  diherence  is  that  here  the  coefficients  of 
pressure  increase  slightly  as  distance  towards  the  tail  of  the  store  increases,  whereas 
in  the  shear  conhguration  they  showed  a  slightly  decreasing  trend.  Table  13  sum- 


marizes  the  average  pressure  coefficient  differential,  the  center  of  pressure,  and  the 
store’s  center  of  mass.  The  pitch-up  mechanism  here  is  similar  to  that  of  the  shear 
configuration,  although  here  the  pressure  differential  is  greater,  resulting  in  a  higher 
pitch- up  moment. 
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Figure  86:  Pressure  coefficient  on  top  and  bottom  of  store  in  forward  shear  config¬ 
uration 


Table  13:  Average  pressure  differential,  center  of  pressure,  and  center  of  mass  for 
the  store  in  forward  shear  conhguration 


Average  Pressnre  Diherential 

Cp^top  C phQtioYn 

Center  of  Pressnre 
x/L 

Center  of  Mass 
x/L 

-0.136 

-0.055 

0.459 

The  mean  and  standard  deviation  of  the  longitudinal  force  and  pitch  moment 
between  the  two  shear  conhgurations  are  compared  in  Table  14.  The  forces  on  the 
store  are  about  equivalent  in  both  cases,  although  the  forward  shear  conhguration 
experiences  a  greater  pitch  up  moment. 


Table  14:  Comparison  of  store  forces  and  moments  between  shear  conhgurations 


Shear  confignration 

Forward  shear  conhgnration 

Mean 

Standard  Deviation 

Mean 

Standard  Deviation 

Longitndinal 

force 

11.19  N 

7.67  N 

11.64  N 

9.01  N 

Vertical 

force 

-0.37  N 

2.52  N 

-0.38  N 

2.68  N 

Pitch 

moment 

0.73  N-m 

2.87  N-m 

1.14  N-m 

3.17  N-m 

The  longitudinal  and  vertical  force  spectra  are  plotted  alongside  the  K16  trans¬ 
ducer  spectrum  in  Figs.  87  and  88.  The  force  spectra  are  plotted  on  a  logarithmic 
scale  in  order  to  match  the  scaling  with  the  cavity  ceiling  SPL.  The  force  spectra  show 
a  strong  correlation  to  the  pressnre  huctnations  on  the  cavity  ceiling.  The  solntions 
here  show  abont  the  same  correlation  to  the  ceiling  hnctnations  as  were  observed  in 
the  shear  conhguration.  The  hrst  fonr  modes  as  observed  on  the  cavity  ceiling  are 
apparent  in  the  store  forces.  The  relative  amplitndes  of  the  peaks  from  the  broad¬ 
band  noise  levels  also  correlate.  The  pitch  and  roll  moment  spectra  are  also  plotted 
alongside  the  K16  transducer  spectrum  in  Figs.  89  and  90.  As  with  the  force  spectra, 
the  moment  spectra  are  plotted  on  a  logarithmic  scale.  Again,  note  that  the  pitch 
and  roll  moments  are  plotted  on  diherent  scales  and  diher  by  about  three  orders  of 
magnitude.  As  was  observed  in  the  carriage  and  shear  conhgnrations,  the  pitch  mo¬ 
ment  spectrnm  shows  a  very  strong  correlation  to  the  K16  SPL  spectrnm.  The  modal 
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frequencies  and  the  peak  relative  amplitudes  match  very  well  with  the  K16  SPL.  The 
roll  moment  spectrum  does  not  show  any  correlation  to  the  K16  spectrum.  None  of 
the  cavity  modes  are  present  in  the  roll  moment  data. 


Cavity  Ceiling  Pressure  and  Store  Longitudinal  Force  Spectra 


o 

CL 


Figure  87:  Store  longitudinal  force  frequency  spectrum  (right  axis)  and  K16  trans¬ 
ducer  SPL  spectrum  (left  axis)  for  forward  shear  conhguration 
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Cavity  Ceiling  Pressure  and  Store  Vertical  Force  Spectra 


Figure  88:  Store  vertical  force  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  forward  shear  conhguration 


Cavity  Ceiling  Pressure  and  Store  Pitch  Moment  Spectra 


Figure  89:  Store  pitch  moment  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  forward  shear  conhguration 
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Cavity  Ceiling  Pressure  and  Store  Roll  Moment  Spectra 
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Figure  90:  Store  roll  moment  frequency  spectrum  (right  axis)  and  K16  transducer 
SPL  spectrum  (left  axis)  for  forward  shear  conhguration 


The  SPL  spectra  on  the  surface  of  the  store  are  presented  in  Figs.  91  through 
93.  Each  of  the  hgures  presents  the  SPL  spectra  at  a  given  longitudinal  position  on 
the  store.  The  pressure  spectra  at  each  longitudinal  position  correlate  well  with  each 
other.  At  x/L  =  0.25,  all  four  modes  are  clearly  present  with  the  fourth  mode  being  the 
dominant  frequency.  At  x/L  =  0.50,  all  four  modes  are  still  clearly  evident,  although 
here  the  second,  third,  and  fourth  modes  are  roughly  equivalent  in  amplitude.  At 
the  x/L  =  0.75  location,  the  hrst,  third,  and  fourth  modes  are  present,  although  the 
third  mode  is  greatly  diminished  here  and  the  second  mode  is  the  dominant  mode. 
The  hrst  mode  here  is  not  observable.  It  is  also  of  note  that  at  all  three  positions, 
the  bottom  and  side  pressure  data  correspond  well  in  overall  sound  levels,  however 
only  at  x/L  =  0.75  does  the  SPL  at  the  top  of  the  store  corresponds  with  the  other 
two  locations.  Here,  though,  all  three  positions  correspond  very  well  with  each  other. 
Note  that  the  carriage  and  shear  conhgurations  showed  similar  trends  in  the  pressure 
loading  across  the  store.  These  same  trends  are  not  observed  here.  This  observation 


93 


suggests  that  the  turbulent  characteristics  of  the  flow  vary  mainly  with  longitudinal, 
rather  than  vertical,  position  in  the  cavity. 


Average  Pressure  at  X/L  =  0.25 


Figure  91:  SPL  spectra  on  store  body  at  x/L  =  0.25  for  forward  shear  conhguration 
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Average  Pressure  at  X/L  =  0.50 


Figure  92:  SPL  spectra  on  store  body  at  x/L  =  0.50  for  forward  shear  configuration 


Average  Pressure  at  X/L  =  0.75 


Figure  93:  SPL  spectra  on  store  body  at  x/L  =  0.75  for  forward  shear  conhguration 
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The  same  vibrational  modes  which  are  observed  in  the  empty  cavity  are  also 
present  for  the  cavity  with  a  store  located  upstream  in  the  shear  layer.  Although 
previous  experiments  have  shown  that  a  store  located  close  to  the  cavity  leading 
edge  can  act  as  a  suppressor,  such  an  effect  was  not  observed  here.  In  fact,  the 
response  at  higher  frequency  modes  was  observed  to  increase  in  this  configuration. 
The  suppression  effect  previously  observed  is  likely  reliant  on  the  store  overall  shape, 
size,  and  position.  A  comparison  of  the  cavity  surface  pressure  fluctuations  and  the 
store  force  and  moment  loadings  shows  correlation  between  the  longitudinal  force, 
vertical  force,  and  pitch  moment  with  the  cavity  vibrational  modes.  The  dominant 
moment  on  the  store  in  the  forward  shear  location  is  a  pitch- up  moment.  The  pressure 
fluctuations  on  the  store  in  the  cavity  vary  by  longitudinal  position  and  correlate 
between  circumferential  position  at  a  given  longitudinal  position.  The  cavity  response 
and  store  force/moment  loading  for  this  configuration  are  consistent  with  what  was 
observed  for  the  previous  two  configurations. 
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V.  Conclusions 


A  computational  investigation  was  completed  in  order  to  investigate  the  relationship 
between  the  acoustic  modes  of  a  cavity  in  grazing  flow  and  the  relationship  of  these 
modes  to  the  force  and  moment  loading  of  a  store  located  in  the  cavity.  A  cavity  of 
L/D  =  4.5  with  freestream  flow  at  M  =  0.95  was  investigated.  Three  conhgurations 
with  a  store  present  were  investigated  and  compared  to  an  empty  cavity:  one  conhgu- 
ration  with  the  store  located  in  the  carriage  position  and  two  conhgurations  with  the 
store  located  in  the  shear  layer.  Several  conclusions  can  be  drawn  from  the  results  of 
this  study. 

First  and  foremost,  it  was  demonstrated  that  the  force  and  moment  loading 
on  the  store  correlate  to  the  pressure  huctuations  on  the  cavity  surface.  Specihcally, 
the  longitudinal  force,  vertical  force,  and  pitch  moment  exhibited  huctuations  at  the 
same  frequencies  as  the  cavity  resonant  modes  and  showed  peaks  at  similar  noise  lev¬ 
els  relative  to  the  broadband  noise  as  the  cavity  ceiling  and  wall  pressure  huctuations. 
The  other  forces  and  moments  were  essentially  zero-valued  and  showed  random  huc¬ 
tuations.  This  correlation  holds  for  the  store  regardless  of  position  within  the  cavity. 
The  turbulent  huctuations  observed  throughout  the  cavity  are  consistent.  A  variation 
in  turbulent  huctuations  was  not  observed  with  a  change  in  vertical  position  when 
longitudinal  position  was  kept  constant,  and  the  same  huctuations  were  observed  in 
the  shear  layer  as  well  as  within  the  cavity  itself.  Essentially,  the  store  integrated 
unsteady  loads  correlate  to  the  acoustic  modes  of  the  cavity. 

For  the  cases  observed  here,  the  presence  of  the  store  in  the  cavity  did  not 
fundamentally  alter  the  acoustic  response  that  was  observed  in  the  cavity.  In  the 
three  conhgurations  investigated  here,  the  vortex  shedding  process  and  the  cavity 
huid  response  were  largely  unimpeded  by  the  presence  of  the  store.  This  observation 
further  aids  the  store  separation  engineer  since  it  implies  that  the  characterization 
of  an  empty  cavity  can  be  used  for  the  analysis  of  the  separation  of  entire  classes  of 
stores  from  that  cavity.  There  is  some  evidence,  though,  that  the  presence  of  the  store 
may  have  dehected  some  of  the  how  down  and  out  of  the  cavity,  which  ahected  the 
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low-frequency  response  at  the  aft  wall.  A  suppression  of  cavity  vibrations  due  to  the 
presence  of  a  store  in  the  shear  layer  was  not  observed  here.  Even  in  the  forward  shear 
conhguration,  the  overall  noise  levels  and  the  peak  amplitudes  were  not  diminished. 
It  is  expected  that  this  effect  depends  greatly  on  specihc  store  geometry  and  position. 

Previous  studies  measured  pressure  data  only  along  the  cavity  centerline  with 
the  intent  of  measuring  the  longitudinal  modes  of  the  cavity.  Here,  data  was  taken  on 
the  sidewalls  as  well.  These  data  did  not  provide  any  evidence  of  any  signihcant  lateral 
modes  for  this  geometry.  The  frequency  spectra  on  the  sidewalls  showed  the  same 
acoustic  modes  as  those  on  the  cavity  ceiling  and  aft  wall.  Additionally,  the  frequency 
content  on  the  sidewall  locations  correspond  well  to  the  corresponding  longitudinal 
locations  on  the  cavity  ceiling  and  to  the  pressure  fluctuations  along  the  store.  Higher 
overall  sound  levels  and  more  low  frequency  oscillations  are  observed  at  downstream 
locations  compared  to  upstream  locations.  The  observed  modes  are  dependent  on 
cavity  geometry,  but  for  deep  cavities  with  the  characteristic  length  in  the  direction 
of  the  flow,  the  longitudinal  modes  dominate  the  acoustic  response  in  the  cavity. 

Finally,  it  was  shown  that  the  overall  force  and  moment  loading  on  a  store  in  a 
cavity  are  a  total  cavity  phenomenon.  In  other  words,  the  same  forces  and  moments 
which  are  signihcant  in  the  carriage  position  are  also  signihcant  in  the  shear  layer. 
Most  signihcantly,  it  was  observed  here  that  the  released  store  is  subjected  to  a  pitch- 
up  moment  from  the  instant  it  is  released  from  carriage.  This  moment  loading  is  a 
result  of  the  mean  how  in  the  cavity;  the  pitch-up  moment  is  not  merely  a  result  of 
the  store  penetrating  the  shear  layer.  However,  since  the  mean  how  throughout  the 
cavity  is  not  uniform,  the  mechanisms  which  create  the  force  and  moment  loading 
are  not  consistent  throughout  the  cavity.  There  is  evidence  to  suggest  that  the  center 
of  pressure  of  the  store  shifts  forward  as  it  moves  downward  through  the  cavity,  and 
the  store  transitions  from  a  stable  conhguration  to  an  unstable  one  as  a  result.  The 
behavior  of  a  store  as  it  is  released  from  a  cavity  may  be  better  explained  not  by  a 
change  in  the  overall  force  and  moment  loading  on  the  store  as  it  enters  the  shear 


layer,  but  rather  a  change  in  response  to  the  store  loading  brought  on  by  a  change  in 
aerodynamic  stability. 

5. 1  Future  Work 

The  results  and  conclusions  presented  here  open  the  door  for  future  areas  of 
research. 

The  conhgurations  investigated  here  are  representative  of  only  one  cavity  geom¬ 
etry  with  one  store  loadout.  Additional  studies  are  necessary  to  fully  characterize  the 
many  different  conhgurations  which  exist.  A  store  which  is  relatively  larger  compared 
to  the  cavity  might  create  a  situation  where  the  store’s  presence  affects  the  cavity’s 
native  response.  Asymmetric  weapon  loadouts  could  also  signihcantly  affect  the  how 
within  the  cavity.  It  was  also  observed  that  the  turbulence  levels  and  frequencies  vary 
within  the  cavity  based  upon  longitudinal  position.  A  further  investigation  into  the 
cavity  response  and  store  force  and  moment  loading  for  diherent  cavity  geometry /s- 
tore  combinations  is  necessary  to  fully  understand  this  problem.  An  in-depth  study 
investigating  the  ehect  of  varying  store  sizes,  shapes,  and  position  within  the  cavity 
is  warranted. 

A  more  robust  study  of  store  stability  within  the  cavity  needs  to  be  conducted. 
The  observations  from  this  research  indicate  that  the  how  conditions  within  a  cavity 
can  cause  a  store  which  is  initially  stable  to  become  unstable  as  it  exits  the  cavity. 
A  better  understanding  of  this  transition  is  necessary  in  order  to  be  able  to  properly 
predict  its  ehect  on  a  store’s  trajectory.  A  store  stability  analysis  can  be  accom¬ 
plished  through  a  combination  of  static  solutions  where  the  position  of  the  store  is 
varied  as  well  as  with  dynamic  cases  where  a  trajectory  is  simulated  and  the  pressure 
distributions  and  center  of  pressure  are  calculated  throughout  the  trajectory.  Fully 
understanding  this  phenomenon  is  crucial  in  order  to  develop  techniques  to  maintain 
safe  and  acceptable  separation  of  cavity-released  stores. 
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It  has  been  demonstrated  that  the  force  and  moment  loading  on  a  store  in  a 
cavity  correlate  to  the  acoustic  modes  of  the  cavity.  The  effects  of  this  relationship 
on  a  released  store’s  separation  characteristics  need  to  be  examined.  This  would 
be  done  by  investigating  the  change  in  a  store’s  trajectory  as  various  parameters  are 
altered,  including  time  of  release  and  store  mass  and  inertial  properties.  As  the  effects 
of  the  store  loading/acoustic  mode  relationship  are  better  understood,  a  predictive 
capability  can  be  better  developed. 

Investigation  into  realistic  aircraft  cavity  geometries  need  to  be  conducted.  Al¬ 
though  much  can  be  learned  from  generalized  geometries  about  the  fundamental  pro¬ 
cesses  at  work,  geometries  which  are  representative  for  operational  aircraft  need  to 
be  studied  in  order  for  this  knowledge  to  be  transferable  to  real-world  applications. 
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Appendix  A.  OVERFLOW  Namelist  Parameters 

This  appendix  tabulates  the  signihcant  OVERFLOW  namelist  parameters  that  were 
used  for  this  research.  These  are  provided  as  a  reference  so  that  the  solutions  here 
can  be  easily  recreated  for  future  research.  The  inputs  are  divided  into  sections  that 
mimic  the  same  groupings  of  variables  as  the  OVERFLOW  2.1  User’s  Guide  [18].  For 
any  variable  not  listed  in  Tables  15  through  23,  the  OVERFLOW  2.1  default  value  is 
used. 

Table  15:  Summary  of  Global  OVERFLOW  (SGLOBAL)  inputs 


Parameter 

Value 

NSTEPS 

10000 

RESTRT 

.F. 

NSAVE 

200 

ISTART.QAVG 

1000 

NQT 

205 

NQG 

0 

MULTIG 

.F. 

FMG 

.T. 

FMGGYG 

150,  150 

NITNWT 

5 

FSONWT 

2 

DTPHYS 

1.93 

Table  16:  Summary  of  Global  OVERFLOW-D  (SOMIGLB)  inputs 


Parameter 

Value 

DYNMGS 

.F. 

I6DOF 

2 

IRUN 

0 

NADAPT 

0 

LFRINGE 

3 

Table  17:  Summary  of  OVERFLOW  DGF  (SDGFGLB)  inputs 


Parameter 

Value 

DQUAL 

0.2 

MORFAN 

1 

NORFAN 

13 
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Table  18:  Summary  of  load  balance  (SGROUPS)  inputs 


Parameter 

Value 

WGHTNB 

2.0 

IGSIZE 

500000 

MAXNB 

0 

MAXGRD 

0 

Table  19:  Summary  of  flow  parameter  (SFLOINP)  inputs 


Parameter 

Value 

ALPHA 

0.0 

BETA 

0.0 

FSMAGH 

0.95 

GAMINF 

1.4 

REY 

20833.3 

PR 

0.72 

PRT 

0.9 

TINF 

518.67 

Table  20:  Summary  of  numerical  method  selection  (SMETPRM)  inputs 


Parameter 

Value 

ILHS 

6 

IRHS 

5 

ILIMIT 

4 

Table  21:  Summary  of  time  accuracy  (STIMACU)  inputs 


Parameter 

Value 

ITIME 

0 

DT 

0.0 

TFOSO 

1.0 

Table  22:  Summary  of  smoothing  parameter  (SSMOACU)  inputs 


Parameter 

Value 

DELTA 

ESO 

1.1 

5 
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Table  23: 


Summary  of  viscous  and  turbulence  modeling  (SVISINP)  inputs 


Parameter 

Value 

Vise 

.T. 

IDES 

2 

IRC 

0 

WALLFUN 

.T. 
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Appendix  B.  Summary  of  Computational  Expense  of  Solutions 

This  appendix  summarizes  the  computational  expense  for  the  OVERFLOW  2.1  so¬ 
lutions  conducted  during  this  research.  In  cases  where  solutions  were  computed  on 
differing  computer  systems,  the  system  used  is  also  noted.  Additional  commentary  is 
also  provided. 

An  overview  of  the  computational  expense  of  the  main  solutions  are  presented 
in  Table  24.  These  solutions  were  all  computed  on  the  AFIT  local  cluster. 


Table  24:  Computational  Cost  of  Cavity  Simulations 


Number 
of  cells 

Cores  Used 

Wall  Time 
(hrs) 

CPU  Time 
(cpu  ■  hrs) 

Empty  cavity 

5,172,360 

50 

53.5 

2677 

Carriage  configuration 

9,357,596 

28 

156 

4368 

Shear  configuration 

9,357,596 

30 

165 

4950 

Forward  shear  configuration 

9,357,596 

50 

96.1 

4804 

The  computational  expense  for  each  of  the  solutions  computed  for  the  Newton 
sub  iteration  study  are  displayed  in  Table  25.  The  percent  change  is  calculated  using 
zero  subiterations  as  the  baseline.  There  is  a  slight  decrease  in  cost  when  using 
one  subiteration  compared  to  zero  subiterations.  Using  Eve  subiterations  increases 
the  computational  cost  by  almost  5  times  and  using  10  subiterations  increases  the 
computational  cost  by  almost  10  times  from  the  baseline.  This  increased  expense 
makes  it  critical  to  use  the  fewest  number  of  subiterations  necessary  for  an  accurate 
answer  in  order  to  minimize  the  computational  cost  of  computing  the  solutions. 


Table  25:  Comparison  of  Computational  Cost  for  Varying  Number  of  Subiterations 


Cores  Used 

Wall  Time 
(hrs) 

CPU  Time 
(cpu  ■  hrs) 

%  Change 

0  subiterations 

30 

14.4 

431 

- 

1  subiteration 

50 

7.0 

351 

-18.6% 

5  subiterations 

30 

78.4 

2352 

446% 

10  subiterations 

50 

92.7 

4637 

976% 

The  computational  expense  for  each  of  the  solutions  computed  for  the  grid 
resolution  study  are  displayed  in  Table  26.  The  percent  change  is  calculated  using 
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the  medium  grid  as  the  baseline.  The  medium  grid  actually  took  more  computational 
time  to  run  than  the  hue  grid,  despite  having  fewer  grid  points.  The  additional  run 
time  is  likely  due  to  the  fact  that  the  medium  grid  solution  was  run  on  the  AFIT  local 
cluster  whereas  the  coarse  and  hne  grid  solutions  were  run  on  the  Air  Force  Research 
Laboratory  Department  of  Defense  Supercomputing  Resource  Center  (AFRL  DSRC) 
system,  which  means  the  comparison  between  these  solutions  is  not  direct.  The  AFRL 
DSRC  system  is  faster  than  the  AFIT  cluster,  which  is  why  the  fine  grid  solution  was 
able  to  be  completed  in  less  computational  time  than  the  medinm  grid  solution  even 
thongh  it  had  more  points  overall. 


Table  26:  Comparison  of  Compntational  Cost  for  Grid  Resolution  Stndy 


Number 
of  cells 

Cores  Used 

Wall  Time 
(hrs) 

CPU  Time 
(cpu  •  hrs) 

%  Change 

System 

Coarse  grid 

4,650,360 

50 

35.4 

1771 

-33.8% 

DSRC 

Medium  grid 

5,172,360 

50 

53.5 

2677 

- 

AFIT 

Fine  grid 

6,476,760 

50 

50.2 

2508 

-6.3% 

DSRC 

The  computational  expense  for  each  of  the  solutions  computed  for  the  grid 
resolution  study  are  displayed  in  Table  27.  The  percent  change  is  calculated  using 
dt  =  1.6  X  10“^  s  as  the  baseline.  As  with  the  grid  sensitivity  solution,  a  direct 
comparison  between  the  cost  of  these  solntions  cannot  be  made.  The  smaller  and 
larger  time  steps  were  rnn  on  the  AFRL  DSRC,  whereas  the  medinm  time  step  was 
run  on  the  AFIT  cluster.  A  comparison  between  the  small  and  large  time  step  shows, 
however,  that  the  choice  of  time  step  does  not  affect  the  compntational  cost  of  the 
solntion. 


Table  27:  Comparison  of  Computational  Cost  for  Time  Step  Sensitivity  Stndy 


dt  (s) 

Cores  Used 

Wall  Time 
(hrs) 

CPU  Time 
(cpu  ■  hrs) 

%  Change 

System 

8.0  X  10“® 

50 

38.9 

1943 

-27.5% 

DSRC 

1.6  X  10“^ 

50 

53.5 

2677 

- 

AFIT 

4.0  X  10-® 

50 

39.0 

1951 

-27.1% 

DSRC 
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